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ABSTRACT 

Thus far, judging the fate of a massive star (either a neutron star (NS) or a black hole) solely by its structure 
prior to core collapse has been ambiguous. Our work and previous attempts find a non-monotonic variation of 
successful and failed supernovae with zero-age main-sequence mass, for which no single stmctural parameter 
can serve as a good predictive measure. However, we identify two parameters computed from the pre-collapse 
stmcture of the progenitor, which in combination allow for a clear separation of exploding and non-exploding 
cases with only few exceptions (~l-2.5%) in our set of 621 investigated stellar models. One parameter is 
M4, defining the normalized enclosed mass for a dimensionless entropy per nucleon of s = 4, and the other is 
yu4 = (d;»/M Q )/(dr/1000km)| s= 4, being the normalized mass-derivative at this location. The two parameters JI 4 
and M 4 IA 4 can be directly linked to the mass-infall rate, M, of the collapsing star and the electron-type neutrino 
luminosity of the accreting proto-NS, L Ve oc M ns M, which play a crucial role in the “critical luminosity” 
concept for the theoretical description of neutrino-driven explosions as runaway phenomenon of the stalled 
accretion shock. All models were evolved employing the approach of Ugliano et al. for simulating neutrino- 
driven explosions in spherical symmetry. The neutrino emission of the accretion layer is approximated by a 
gray transport solver, while the uncertain neutrino emission of the 1.1 M 0 proto-NS core is parametrized by an 
analytic model. The free parameters connected to the core-boundary prescription are calibrated to reproduce 
the observables of Supernova 1987A for five different progenitor models. 

Subject headings: supernovae: general — stars: massive — hydrodynamics — neutrinos 


1. INTRODUCTION 

Presupernova stars in the mass range above ~9 M 0 exhibit 
large variations of their structure with respect to, e.g., their 
Fe-core and O-core masses, their binding energies and their 
density or entropy profiles above the Fe-core (Woosley et al. 
2002). These properties vary non-monotonically with the 
zero-age main-sequence (ZAMS) mass and can differ consid¬ 
erably even between progenitors with onl y a s mall difference 
of their ZAMS masses (fSukhbold & Woosle y|2014| i. 


Correspondingly, [Ugliano et al. ~|2012) found that the 
properties of neutrino-driven supernovae (SNe) like explo¬ 
sion energy, nickel mass, and remnant mass change non- 
monotonically with the ZAMS mass. In particular, for the in¬ 
vestigated grid of 101 solar-metallicity progenitors binned in 
0.2 M 0 steps (Woosley et al. |2002| ), they found islands of non¬ 
exploding, black hole (BH) forming cases down to 15 M 0 , al¬ 
ternating with mass-intervals of exploding progenitors. In a 
few cases individual neighboring progenitors showed oppo¬ 
site behavior. 

Ugliano et al. (2012[> used a simple, parametric model for 


the contracting proto-neutron star (PNS) as a neutrino source 
to trigger neutrino-driven explosions in spherically symmet¬ 
ric (ID) hydrodynamic simulations, but their basic findings 
were confirmed by other groups working with semi-analytic 
descriptions in ID (Pejcha & Thompson |2015 1 and approx¬ 
imate neutrino transport in two- and three-dimensional (2D, 
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3D) hydrod ynamic models (Naka mura et al. 2015 Horiuchi 
|et al.||2014|). While |0’Connor & Ott|(|201 1 1 suggested that 
BH formation requires a compactness (normalized enclosed 
mass-radius ratio) of £2.5 > 0.45 with 


(tm = 


M/M q 


R(M)/ 1000 km ’ 


( 1 ) 


Ugliano et al. 
explosions or 


( |2012| > obtained only explosions for £2.5 < 0.15, 
BH formation for 0.15 < £2.5 <0.35, and only 
BH formation for £2.5 > 0.35, which implies a larger frac- 
tion of BH forma tion cases for solar-metallicity stars. Hori-| 
uchi et al. (j2014f pointed out that a critical compactness of 
£2.5 ^ 0.2 for failed explosions is compatible with a lack 
of red supergiant Type-IIP SN progenitors above ~16M 0 
( Smartt et al.|2009 1 and with a significant excess o f the star- 
format ion r ate compared to the observ ed SN rate ( |Horiuchi1 
et al. 201 H >. Pejcha & Thompson ( 2015| > showed that their pa¬ 
rameterization “c ase (a) ”, which yields results similar to those 
of Ugliano et al. (j2012jl, is close to being optimally compati¬ 
ble with a combination of several observational constraints. 

How can the non-monotonicities of the explodability be un¬ 
derstood in terms of the pre-supernova properties and in the 
context of the physics of the neutrino-driven mechanism? Are 
there characteristic parameters of the pre-supernova star that 
decide better about success or failure of the explosion than a 
single value of the compactness or other, similarly useful pa¬ 
rameters like the iron-core mass or the binding energy outside 
of the iron core? While all these measures reflect trends like 
an enhanced tendency of BH formation for high compactness, 
large iron-core mass or high exterior binding energy, there are 
still many outliers that do not obey the correlations. For exam¬ 
ple, a suitably chosen mass M of the compactness £m allows 


to correctly predict explosions in <90% of the cases (Pejcha 
& Thompson|2015), but the best choice of M is merely empir- 
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ical and the physical justification of £m as a good diagnostics 
is unclear. 

Here we propose a two-parameter criterion that separates 
successful explosions from failures with very high reliability. 
While two compactness values, e.g. t;\ 5 and £ 2.5 > or the iron- 
core mass and the mean entropy in some suitable mass range 
begin to show such a disentanglement, we demonstrate that 
the normalized mass inside a dimensionless entropy per nu¬ 
cleon of s = 4, 

M 4 = m(s = 4) /M 0 , (2) 


and the mass derivative at this location, 

dm /M q 


P4 = 


dr /1000 km 


(3) 


s=4 


both determined from the pre-supernova profiles, allow to pre¬ 
dict the explosion behavior successfully in > 97% of all cases 
and have a direct connection with the theoretical basis of the 
neutrino-driven mechanism. 

We briefly describe our numerical approach in Sect. [2] in¬ 
cluding a detailed discussion of our modeling methodology 
in comparison to other approaches in the recent literature, 
present our results in Sect. [5] and conclude in Sect. [4] 

2. NUMERICAL SETUP AND PROGENITOR MODELS 
2.1. Modeling approach 

Our basic modeling approach follows Ugliano et al. ( |2012| 
with a number of improvements. To trigger neutrino-driven 
explosions in spherically symmetric (ID) hydrodynamic sim¬ 
ulations, we use a schematic model of the high-density core 
of the PNS as neutrino source (for details, see Ugliano et al. 
2012). This analytic description is applied to the innermost 
1.1 M g , which are excised from the computational domain, 
and it yields time-dependent neutrino luminosities that are im¬ 
posed as boundary values at the contracting, Lagrangian inner 
grid boundary. On the numerical grid, where neutrino opti¬ 
cal depths increase from initially ~10 to finally several 1000, 
neutrino transport is ap proximated by the gray treatment de¬ 
scribed in |Scheck et al.| ( [2006| > and |Arcones et al.| ( |2007 1. This 
allows us to account for the progenitor-dependent variations 
of the accretion luminosity. 

Our approach replaces still uncertain physics connected to 
the equation of state (EoS) and neutrino opacities at high den¬ 
sities by a simple, computationally efficient PNS core model. 
The associated free parameters are calibrated by reproducing 
observational properties of SN 1987A. We emphasize that the 
neutrino emission is sensitive to the time and progenitor de¬ 
pendent mass accretion rate. Not only the accretion luminos¬ 
ity increases for progenitors with higher mass accretion rate 
of the PNS, but also the neutrino loss of the inner core rises 
with the accreted mass because of compressional work of the 
accretion layer on the core. Such dependences are accounted 
for in our modeling of NS core and accretion. 

Our numerical realization improves the treatment by 




Fig. 1.— Electron fraction as function of density, Y e {p), for modeling the 
deleptonization of the collapsing stellar iron core during the infall p has e unti l 
core bounce according to the approximate treatment of |Liebendorfer|(2005f . 

in ejecta with Y e < 0.49 and thus tracks the ejection of mat¬ 
ter with neutron excess, when detailed nucleosynthesis cal¬ 
culations predict little production of 56 Ni (Thielemann et al. 
179961 . 

The network is consistently coupled to the hydrodynamic 
modeling and allows us to include the contribution from ex¬ 
plosive nuclear burning to the energetics of the SN explo¬ 
sions. The collapse phase until core bounce is modeled with 
the deleptonization scheme proposed by |Liebendorfer[|2005[ ), 
using the Y e (p ) trajectory of Fig. [T] for the evolution of the 
electron fraction Y e as function of density p (jMiiller||2013 
private communication). This yields good overall agreement 
with full neutrino transport results and allows for a very effi¬ 
cient computation of large sets of post-bounce models. 

2.2. Progenitor models 

We perform collapse and explosion simulations for large 
progenitor sets of different metallicities, namely: the zero- 
metallicity z2002 set (30 models with ZAMS masses of 11.0- 
40.0 M 0 ), low-metallicity (10~ 4 solar) u2002 series (247 mod¬ 
els, 11.0-75.0 Mq), and the solar-metallicity s2002 series 
(101 models, 10.8-75.0 M 0 ) of Woosley et al.j (j2002| plus a 
10.0 M q progenitor (Woosley 2007 ] pri vate communication) 
and a 10.2 M 0 progenitor ( jHeger||2003| private communica¬ 
tion); the solar-metallicity s20fr(l51 models, 15.0-30.0 M 0 ) 
and sh2014 series (15 models, 30.0-60.0 M 0 , no mass loss) of 
Sukhbold & Woosley| (|2014|, supplemented by additional 36 
models with 9.0-14.9 M 0 ; the solar-metallicity s2007 series 
(32 models, 12.O-12O.OM 0 ) of |Woosley et al. ( 2007 1 ; and the 
n2006 series (8 models, 13.0-50.0 M 0 ; Nomoto et a ,|20061 l. 

For the core-model parameter calibration we choose five 
different progenitors, namely the (red supergiant) model sl9.8 
of the s2002 series as in Ugliano et al. (2012 1 , and four 


blue supergiant pre-supernova models of SN 1987A: wl5.0 


tracer nucleus of iron-group material) at lower temperatures 
([Muller|[l986|. The tracer nucleus is assumed to be formed 


(ZAMS mass of 15 M 0 ; Woosley et al.|1988j >. wl8.0 (18 M 0 , 
evolved with rotation; | Woos ley et al.||2007j >, w20.0 (20 M 0 ; 
Woosley et al.| |1997|l, and n20.0 (20 M 0 ; |Shigeyama & 
Nomoto|1990p . Compactness values and explosion and rem- 
nant parameters of these models are listed in Table [I] 

The calibration aims at producing the explosion energy and 
ejected 56 Ni mass of SN 1987A compatible with observations, 
for which the best values are E exp = (1.50 + 0.12) x 10 51 erg 
(Utrobin 2005 i, E exp ~ 1.3 x 10 51 erg (Utrobin & Chugail 




























































































Explosion criterion of massive stars 


3 


TABLE 1 

Calibration models with explosion and remnant properties 


Calibration 

Model 

fl.5 a 

fl.75 a 

?2.0 a 

^2.5 3 

t b 
‘exp 

[ms] 

Eex P c 

[B] 

M ej d 

[M 0 ] 

^exp / 

[B /M q \ 

[M 0 ] 

^tracer ^ 

[M Q ] 

M ns e 

[M 0 ] 

^wind 

[M Q \ 

Mft 1 

[M G ] 

tv, 90 J 

[s] 

Z7 k 

c v,tot 

[100 B] 

sl9.8 (2002) 

1.03 

0.35 

0.22 

0.14 

750 

1.30 

12.98 

0.100 

0.072 

0.034 

1.55 

0.096 

0.00298 

4.27 

3.68 

wl5.0 

0.34 

0.09 

0.03 

0.01 

580 

1.41 

13.70 

0.103 

0.045 

0.046 

1.32 

0.088 

0.00018 

5.18 

2.81 

wl8.0 

0.76 

0.26 

0.16 

0.10 

730 

1.25 

15.42 

0.081 

0.056 

0.036 

1.48 

0.081 

0.00310 

4.16 

3.32 

w20.0 

0.98 

0.35 

0.18 

0.06 

620 

1.24 

17.81 

0.070 

0.063 

0.027 

1.56 

0.089 

0.00168 

4.73 

3.61 

n20.0 

0.87 

0.36 

0.19 

0.12 

560 

1.49 

14.84 

0.100 

0.036 

0.052 

1.55 

0.117 

0.00243 

3.97 

3.48 


a Compactness evaluated for a central density of 5 X 10 10 g cm 3 . For wl5.0 the value at core bounce is given because earlier data are not available. 
b Post-bounce time of onset of explosion, when shock expands beyond 500 km (1 B = 1 bethe = 10 51 erg). 
c Final explosion energy, including binding energy of preshock progenitor. 

^ Mass ejected in the explosion. 

e Ejected 56 Ni mass produced by explosive burning with late-time fallback taken into account. 

^ Mass of neutron-rich tracer nucleus ejected in neutrino-driven wind material with neutron excess (fallback is taken into account). 
g Final baryonic neutron-star mass including late-time fallback. 

b Neutrino-driven wind mass measured by mass between gain radius at f e xp and preliminary mass cut before fallback. 

J Fallback mass. 

J Emission time for 90% of the radiated neutrino energy, 
k Total radiated neutrino energy. 


|20TT) , and M Ni = 0.0723-0.0772 M 0 ( [Utrobin et akf20l4) , 
but numb ers reported by othe r authors cover a considerable 
range (cf. |Handy et al.|]2014|for a compilation). The explo¬ 
sion energy that we accept for a SN1987A model in the cal¬ 
ibration process is guided by the ejected 56 Ni mass (which 
fully accounts for short-time and long-time fallback) and a ra¬ 
tio of £exp to ejecta mass in the ballpark of estimates based on 
light-curve analyses (cf. Table [l]for our values). 

Because of the “gentle” acceleration of the SN shock by 
the neutrino-driven mechanism (also in 3D simulations, see 
|Utrobin et al.|2014] >, it is difficult to produce this amount of 
ejected 5b Ni just by shock-induced explosive burning. Ms6 Ni 
in Table [1] mainly measures this component but also contains 
56 Ni from proton-rich neutrino-processed ejecta. However, 
also neutrino-processed ejecta and the neutrino-driven wind 
with a slight neutron excess could contribute significantly to 
the 56 Ni production. The electron fraction Y e of these ejecta 
is set by v e and v e interactions and depends extremely sensi¬ 
tively on the properties (luminosities and spectra) of the emit¬ 
ted neutrinos, which our transport approximation cannot re¬ 
liably predict and which also depend on subtle effects con¬ 
nected to multi-dimensional physics and neutrino opacities. 
For these reasons we consider the 56 Ni as uncertain within 
the limits set by the true 56 Ni yield from our network on the 
low side and, in the maximal case, all tracer material added to 
that. We therefore provide as possible 56 Ni production of our 
models the range Ms6 Ni < < Ms6 Ni + M tracer . 

Different from Ugliano et al. ( 2012) we reduce the com¬ 
pression parameter of the NS core model by a linear relation 
£ ^ (n. 75 ,b (the value of i s measured at core bounce) for 
progenitors with <13.5 M 0 , i.e., we use the function 




£l.75.b 

0.5 


for M < 13.5 M 0 


(4) 


with £ being the value determined from the SN 1987A cali¬ 
bration for a considered progenitor model of that supernova. 
We note that the values of £i.75,b are less than 0.5 for all pro¬ 
genitors below 13.5 M 0 and close to 0.5 for M ~ 13.5 M 0 , for 
which reason Eq. (HI) connects smoothly to the C value applied 
for stars above 133 M 0 according to the SN 1987A calibra¬ 
tion. 

The modification of Eq. ([4]) accounts for the reduced bur¬ 


den of the small mass of the accretion layer of these stars 
with their extremely low compactnesses. Such a modifica¬ 
tion allows us to reproduce the trend to weak explosions ob¬ 
tained in sophisticate d 2D and 3D simulations for low-mass 
iron-core progenitors (Janka et al.|2012||Melson et al.|2015a; 
Muller|2015 1 . We point out that the Crab supernova SN 1054 
is conside red to be connected to the explosi on of a ~1OM 0 
star (e.g., |Nomoto et al.||1982| |Smith|[20l3 1 , and its explo¬ 


sion energy is estim ated to be up to only ~ 1 (P° erg (e.g., Yang 
|& Che valier 2015]). This fact lends support to the results of 
recent, self-consistent ID and mu lti-dimensional supernova 
models of < 10 M 0 stars (e.g.,|Kitaura et al.|2006 ; Fischer et al. 
2010| Melson et al.|2015a] ), whose low explosion energies and 


low nickel production agree with the Crab observations. 

Although as a consequence of our ( reduction the explo¬ 
sion times, fexn, tend to be late for stars in the 10.5-12.5 M 0 
range (Fig. Hb, this behavior also seems to be compatible 
with self-consistent, multi-dimensional simulations of stellar 
explosions in the 11-12M 0 range by |Miiller et al.| (|2013J), 
Muller|( f20151 ), and |Janka et al.| ( [20T2] ), where these stars were 
found to have a long-lasting phase of accretion and simul¬ 
taneous mass outflow after a quite inert onset of the explo¬ 
sion. A more detailed discussion and justification of our mod¬ 
ified treatm ent of low-ZAMS mass cases will be provided in 
Sect. |2.3.4| and can also be found in a follow-up paper by 
[Sukhbold et al.[ ( |20151 >, where the values of all PNS core¬ 
model parameters are tabulated for all calibrations. 

It must be emphasized that in the context of the present 
work the detailed treatment of the low-mass stars is not overly 
important. These stars usually explode fairly easily, indepen¬ 
dent of the treatment of the PNS core model with the original 
or with our revised calibration. Therefore these stars lie far 
away from the boundary curve that separates exploding from 
non-exploding models and whose determination will be our 
main goal in Sect. [3] For this reason exactly the same sepa¬ 
ration line is obtained when the core-model parameter values 
from the SN 1987A calibration are applied to all stars. 

In Table |T] and the rest of our paper, time-dependent struc¬ 
tural parameters of the stars (like compactness values, /./ 4 of 
Eq. [3] the iron-core mass Mp e ) are measured when the stars 
possess a central density of 5 x 10 l() gcm~ 3 , unless otherwise 
stated. This choice of reference density defines a clear stan¬ 
dard for the comparison of stellar profiles of different progen- 
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itors (cf. Appendix A of Buras et al.| |2006 1 . Different from 
the moment of core bounce, which was used in other works, 
our reference density has the advantage to be still close to the 
initial state of the pre-collapse models provided by stellar evo¬ 
lution modeling and therefore to yield values of the structural 
parameters that are more similar to those of the pre-collapse 
progenitor data. Our calibration model wl5.0, however, must 
be treated as an exception. Because pre-collapse profiles of 
this model are not available any more, all structural quantities 
for this case are given (roughly) at core bounce. 

2 . 3 . Methodology and theory of neutrino-driven explosions 
2 . 3 . 1 . Status of “ab initio ” supernova modeling 

“Ab initio”, fully self-consistent simulations of stellar core 
collapse with state-of-the-art treatment of microphysics and 
neutrino transport do not lead to explosions in spherical sym¬ 
metry except for stars with O-Ne-Mg and Fe-cores near the 
low-mass end of SN progenitors (Kitaura et al. |2006 Janka 


|et al. |2008l|2012||Fischer et al.||2010[|Melson et al.| |2()15a). 

2D simulations in the recent past have produced successful 
explosions and underline the fundamental importance of mul¬ 
tidimensional effects, but the true meaning of these results 
with respect to the neutrino-driven mechanism is not finally 
clear and the current situation is diffuse and contradictive. 

On the one hand, some of the 2D explosions set in relatively 


late and might remain on the weak side (e.g., Marek & Janka 
2009] |Suwa||2012| |Suwa et al.||20l4t |Mtiller et al.||2012a|b] 


2013] iMuller & Janka||20141 |Miiller||20151 |Hanke||2014| ), al- 

though such apprehension is speculative because not all sim¬ 
ulations could be continued until the explosion energy had 
saturated (Miiller| |2015 1. On the other hand, the Oak Ridge 
group obtained explosions much earlier after bounce with 
shock evolutions being astonishingly similar for 12, 15, 20, 
and 25 M Q st ars and explosion energies fairly compatib le with 
observations ( |Bruenn et al.|201 3] |2014] >. In contrast, |Dolence| 
et al.| (|2015)> did not find any successes in 2D simulations of 
the same progenitors but they used a different treatment of 
gravity, hydrodynamics, equation of state, neutrino transport, 
and neutrino opacities. The exact reasons for the different 
findings will have to be clarified by detailed tests and com¬ 
parisons. The situation is even more diffuse because current 
3D simulations agree in showing slower explosions compared 


to 2D calcu l ations or even no explosions (e.g., |Hanke et al. 
2012 2013j Tamborra et al.||2014[ Mezzacappa et al .||2015 ' 


Couch|2013[|Couch & 0’Connor|2dl4[ |Takiwaki et al ,|20 1 f 

although some studies have proclaimed the opposite behavior 


( Nordhaus et al.|2010] Burrows| |2013| [Dolence et al.||2013| l. 

So far only a few recent 3D calculations with highly refined 
neutrino treatment have obtained successful shock revival by 
the neutrino-driven mechanism ( |Melson et al.|2015a|b[ |Lentz| 
et al.|20f5) . Interestingly, the 3D simulation of a low-mass 
(9.6 M 0 ) progenitor with detailed neutrino physics, whose ex¬ 
plosion energy approached its saturation level, was found to 


explode more energetically in 3D than in 2D (Melson et al. 


2015a|l. This result is in line with a 2D-3D comparison in 
the 11-12 M g range conducted by |Miiller| (|2015)>. In both 
studies accretion downflows and the re-ejection of neutrino- 
heated matter were observed to be different in 2D and 3D 
because of geometry dependent differences of the Kelvin- 
Helmholtz instability and flow fragmentation. The 3D models 
therefore suggest that explosions in 2D are massively affected 
by the assumption of rotational symmetry around the polar 
grid axis and by an inverse turbulent energy cascade, which 


tends to amplify energy on the largest possible scales (Hanke 
et al.|2012|l and also produces numerical artifacts in the post¬ 
explosion accretion phase of the neutron star ( Miil ler|2015 i. It 
must therefore be suspected that the early onset of explosions 
and the extremely unipolar or bipolar deformations along the 
symmetry axis obtained in many 2D models could be artifacts 
of the imposed symmetry constraints. 

2 . 3 . 2 . Modeling recipes in recent literature 

Before 3D modeling will have become a routine task and re¬ 
sults will have converged, neutrino-driven explosions of large 
sets of progenitor stars can be explored for their observa¬ 
tional implications only by referring to simplified modeling 
approaches. Several different recipes have been introduced 
for this recently. Ugliano et al. (2012 1 used an analytic PNS 
core-cooling model in connection with a neutrino transport 
approximation in ID hydrodynamic explosion simulations (as 
briefly summarized in Sect. [2], thus improving the simpler, 
time-dependent boundary neutrino luminosity prescribed by 
previous users of the simulation code ( Kifonidis et al.||2006[ 
2006] 120081 lArcones et al.||2007| |Ar-| 


Scheck et al. 


2004 


cones & Janka 


2011 1 and the even simpler neutrino light-bulb 
treatment (without any transp ort approximation) applied by 
Janka & Muller | ( fl996] i and |Kifonidis et al-l ( |2003| ). |Q’Connor| 
& Ott (2011 1 resorted to a scaling parameter /h ea t to artifi¬ 


cially enhance the neutrino heating by charged-current pro¬ 
cesses behind the stalled shock in ID hydrodynamic models 
with approximate neutrino treatment. |Nakamura et al.[ ( |2015| 
performed an extensive set of 2D simulations with simpli¬ 
fied neu trino transport despite the grains of salt mentioned 
in Sect. |2.3.1| (see also Horiuchi et al. 2014 for cautioning 
against the 2D results). Pejcha & Thompson] (120151) applied 
a semi-analytic model to determine the onset times of the ex¬ 
plosions, using neutrino luminosities from ID calculations of 
accreting PNSs, and estimated explosion properties by ana¬ 
lytic arguments. [Suwa et al.](j2014|) also performed 2D simu¬ 
lations and suggested analytic approximations for describing 
diffusion and accretion components of the neutrino luminosi¬ 
ties from PNSs and a free-fall treatment for the collapse of 
the overlying stellar layers. Perego et al.] ( [20151 invented a 
method they named “PUSH”, which they applied to trigger 
explosions artificially in their general relativistic, ID hydro- 
dynamic core-collapse and PNS formation modeling with so¬ 
phisticated neutrino transport. PUSH gradually switches on 
and off additional neutrino heating of chosen strength dur¬ 
ing a chosen period of time. This procedure is assumed to 
mimic the effects of multi-D hydrodynamics in the postshock 
region. The extra heating is coupled to the heavy-lepton neu¬ 
trino emission from the PNS. 

All of these recipes contain larger sets of parameters and 
degrees of freedom, which are either varied in exploring dif¬ 
ferent cases (e.g., [Pejcha & Thompso n|20l5| or are adjusted 
by co mparison to more complete models (e.g., |Suwa et al.| 
2014|> or by reprod ucing observational benchmarks like thos e 


2012] Perego et a kpOB] ). 

21 and those in the present 


set by SN 1987A (Ugliano et al. 

The models of Ugliano et al.|(|20l2| and those in the present 
paper assume that the explosion trigger is tightly coupled to 
the physics that reflects the main differences between different 
progenitor stars, namely to the post-bounce accretion history 
of the collapsing stellar core and the corresponding accretion 
luminosities of v e and v e . It will have to be seen whether this 
important aspect of the models remains being supported by 
future developments towards a more complete understanding 
of the physics of the central engine that powers the explosion 
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s2013 14.0 M calib. w18.0 s2013 21.0 M calib. w18.0 s2013 27.0 M calib. w18.0 



Time [s] Time [s] Time [s] 

Fig. 2.— Models sl4.0 (left), s21.0 (middle), and s27.0 (right) of the s2014 progenitor series as exemplary cases of successful explosions with the wl8.0 
calibration. The top panels display as functions of post-bounce time the radius of the outgoing shock (black line), the mass accretion rate measured at 500 km 
(blue line; scale on the right side), and the radii of iron core (orange), M 4 = m(s = 4) (red), M 4 + 0.3 (red dashed) and trajectory of the final mass cut (after 
completion of fallback; purple). The second panels from top show the time evolution of the luminosities of v e , v e and a single species of heavy-lepton neutrinos 
v x as labelled in the plot, measured at 500 km (solid lines) and at the inner grid boundary (dashed lines). The third panels from top show the mean energies of all 
neutrino kinds as radiated at 500 km. The vertical dotted lines indicate the onset time of the explosion as the moment when the outgoing shock passes the radius 
of 500 km. The bottom panels provide the time evolution of the diagnostic energy of the explosion (integrated energy of all postshock zones with positive total 
energy; blue line). Also shown are the kinetic energy (red), gravitational energy (black), and internal energy (orange) as integrals over the whole, final SN ejecta 
between the final mass cut (after fallback) on the one side and the stellar surface on the other. The total (binding) energy (purple) as the sum of these energies 
ultimately converges to the diagnostic energy and both of these energies asymptote to the final explosion energy. While this convergence is essentially reached 
after ~4 s in the case of sl4.0, the expansion of shocked matter in the s21.0 model and thus the energy evolution is slowed down at ~2.7 s by the high densities in 
the stellar core. The s27.0 model becomes gravitationally unbound (i.e., the total binding energy becomes positive) even more slowly because of the very massive 
stellar core. The convergence of total energy and diagnostic energy takes tens of seconds in this case. 
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in the context of the neutrino-driven mechanism. 


2 . 3 . 3 . Motivation of modeling assumptions of this work 


The analytic NS-core model introduced by Ugliano et al. 


(j2012| in combination with their approximate transport solver 
for treating the accretion component of the neutrino luminos¬ 
ity as well as neutrino cooling and heating between PNS and 
shock, is an attempt to realize the tight coupling of accretion 
behavior and explodability in close similarity to wh at is found 
in current 2D simulations (e.g. those of|Marek & Janka|2009' 


Mu ller et al.||2012a|b| |2013 |Miiller & Janka||2014| i. Since 
ID models with elaborate neutrino physics and a fully self- 
consistent calculation of PNS cooling miss the critical condi¬ 
tion for explosions by far, it is not the goal of Ugliano et al. 
( ,2012 1 to closely reproduce the neutrino emission properties 
of such more sophisticated calculations. Rather than that it 
is the goal to approximate the combined effects of neutrino 
heating and multi-dimensional postshock hydrodynamics by 
a simple and computationally efficient neutrino source model, 
which allows for the fast processing of large progenitor sets 
including the long-time evolution of the SN explosion to de¬ 
termine also the shock breakout and fallback evolution. 

Free parameters in the NS core model and the prescribed 
contraction behavior of the inner grid boundary are calibrated 
by matching basic observational features (explosion energy, 
56 Ni yield, total release of neutrino energy) of SN 1987A. 
This is intended to ensure that the overall properties of the 
neutrino-source model are anchored on empirical ground. Of 
course, the setting of the parameter values cannot be un¬ 
ambiguous when only a few elements of a single observed 
SN are used for deriving constraints. However, the approxi¬ 
mate nature of the neutrino source treatment as a whole does 
not require the perfectly accurate description of each individ¬ 
ual model component in order to still contain the essence of 
the physics of the system like important feedback effects be¬ 
tween accretion and outflows and neutrinos, which govern 
the progenitor-dependent variations of explodability and SN 
properties. A reasonable interplay of the different components 
is more relevant than a most sophisticated representation of 
any single aspect of the neutrino source model. 

In detail, the basic features of our ID realization of the 


neutrino-driven mechanism along the lines of Ugliano et al. 
(2012) are the following: 


• The possibility of an explosion is coupled closely to 
the progenitor-dependent strength and evolution of the 
post-bounce accretion. This is achieved not only by tak¬ 
ing into account the accretion luminosity through the 
approximate neutrino transport scheme but also through 
the response of the PNS-core to the presence of a hot 
accretion mantle. The evolution of the latter is explic¬ 
itly followed in our hydrodynamic simulations, which 
track the accumulation of the collapsing stellar matter 
around the inner PNS core. The existence of the man¬ 
tle layer enters the analytic core model in terms of the 
parameter m acc for the mass of this layer and the corre¬ 
sponding accretion rate m acc . 


• The inner 1.1 M 0 core of the PNS is cut out and re¬ 
placed by a contracting inner grid boundary and a corre¬ 
sponding boundary condition in our model. This inner 
core is considered to be the supranuclear high-density 
region of the nascent NS, whose detailed physics is still 
subject to considerable uncertainties. This region is re¬ 
placed by an analytic description, whose parameters T, 


R c (t), and n (see Ugliano et al. 12012b are set to the 
same values for all stars. This makes sense because the 
supranuclear phase is highly incompressible, for which 
reason it can be expected that the volume of the core 
is not largely different during the explosion phase for 
different PNS masses. Moreover, the neutrino diffusion 
time scale out of this core is seconds, which implies that 
its neutrino emission is of secondary importance during 
the shorter post-bounce phase when the explosion de¬ 
velops. Despite its simplicity, our core treatment still 
includes progenitor and accretion dependent variations 
through the mass w acc of the hot accretion mantle of 
the PNS and the mass accretion rate m acc , whose influ- 
en ce on the inner core i s accounted for in Eqs. (1)—(4) 
of Ugliano et al. (2012 1 for describing the energy evo¬ 
lution of the core model. 


The onset of the explosion is considerably delayed (typ¬ 
ically between several 100 ms and about a second) with 
a slow (instead of abrupt) rise of the explosion energy 
during the subsequent shock acceleration phase, when 
an intense neutrino-driven wind ejects matter and de¬ 
livers power to the explosion. Neutrino heating cannot 
deposit the explosion energy impulsively, because the 
ejected matter needs to absorb enough energy from neu¬ 
trinos to be accelerated outwards. The rate of energy 
input to the explosion is therefore limited by the rate at 
which matter can be channeled trough the heating re¬ 
gion. A long-lasting period (hundreds of milliseconds 
to more than a second) of increasing energy is charac¬ 
teristic of neutrino-driven explosions (Fig. [2]) and is ob¬ 
served as gradual growth of the explosion energy also 
in 2 D explosion models, e.g., by |Scheck et al. |2006) 
and |Bruenn et al.| (|2014). To achieve this behavior in 
our ID models the core-neutrino source needs to keep 
up high neutrino luminosities for a more extended pe¬ 
riod of time than found in fully self-consistent SN sim¬ 
ulations, where the rapid decline of the mass-accretion 
rate at the surface of the iron core and at the interface 
of silicon and silicon-enriched oxygen layers leads to a 
strong decrease of the accretion luminosity. The longer 
period of high neutrino emission is compensated by a 
somewhat underestimated early post-bounce neutrino 
luminosity (Fig. [2]i in order to satisfy the energy con¬ 
straints set by the total gravitational binding energy of 
the forming NS. 


• The time scale and duration of the growth of the explo¬ 
sion energy in multi-D models of neutrino-driven SNe 
are connected to an extended period of continued ac¬ 
cretion and simultaneous shock expansion that follows 
after t he revival of the stalled shock (see Ma rek& Janka| 
2009). Persistent accretion thereby ensures the mainte¬ 
nance of a significant accretion luminosity, while par¬ 
tial re-ejection of accreted and neutrino-heated matter 
boosts the explosion energy. In our ID simulations the 
physics of such a two-component flow cannot be accu¬ 
rately accounted for. In order to approximate the conse¬ 
quences of this truly multi-dimensional phase, our ID 
models are constructed with important two properties: 
On the one hand they refer to a high level of the PNS- 
core luminosity for about one second. On the other 
hand they are set up to possess a more extended ac¬ 
cretion phase that precedes the onset of the delayed ex- 
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Fig. 3.— Explosion properties for models of the s2014 series, the supplementary low-mass progenitors with Mzams <15 M©, and models with Mzams > 30 M© 
from the s2007 series, exploded with the wl8.0 calibration. A black vertical line marks the boundary between the two progenitor sets. The panels show the 
final explosion energies, Eexp ( top ; 1 B = 1 bethe= 10 51 erg), times of the onset of the explosion, r e xp (defined as the moment when the shock expands beyond 
500 km; second from top), masses of ejected explosively produced 56 Ni (red bars) and tracer element (orange bars; third from top), baryonic remnant masses with 
fallb ack masses indicated by orange sections of the bars {fourth from top), fallback masses (plotted logarithmically from top), gravitational remnant masses 
(Eq. {TO}: sixth from top), and total energies radiated in neutrinos, £V,tot {bottom). Mass and parameter values of the calibration models are indicated by vertical 
and horizontal blue lines, respectively, with the dashed horizontal blue line in the middle panel giving the 56 Ni mass and the solid horizontal blue line the sum of 
56 Ni and tracer. Non-exploding cases are marked by short vertical black bars in the upper half of each panel. 
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ZAMS Mass [M s ] 


15 20 

ZAMS Mass 


Fig. 4.— Compactness £ 2.5 versus ZAMS mass for the s2002 (left) and s2014 progenitor series (plus Mza ms <15 M Q models) with exploding (red bars) and 
non-exploding (gray bai*s) cases for all calibrations. The s2002 results for the si9.8 calibration agree well with|Ugliano et al.|(2012). Blue vertical and horizontal 
lines indicate the values of the progenitors used for the calibrations. 
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plosion before the intense neutrino-driven wind pumps 
energy into the explosion. The power and mass loss in 
this wind are overestimated compared to sophisticated 
neutrino-cooling simulations of PNSs. However, this 
overestimation of the wind strength has its justification: 
The early wind is supposed to mimic the mass ejection 
that is fed in the multi-dimensional case by the inflow 
and partial re-ejection of matter falling towards the gain 
radius during the episode of simultaneous accretion and 
shock expansion. The enhanced wind mass counterbal¬ 
ances the extra mass accretion by the PNS during the 
long phase before the shock acceleration is launched, 
and this enhanced wind mass is of crucial importance 
to carry the energy of the neutrino-powered blast. 

Figure[2]shows the post-bounce evolution of the stalled SN 
shock, the onset of the explosion, energy evolution, and the 
time evolution of the neutrino emission properties for three 
representative progenitors, namely sl4, s21, and s27 of the 
s2014 series, which explode successfully with the wl8.0 cal¬ 
ibration. The shock stagnation at a radius of approximately 
200 km lasts between ~700 ms and 900 ms and can exhibit 
the well-known oscillatory expansion and contraction phases, 
which signal proximity to the explosion (see e.g., Buras et al. 
2006[ Murphy & Burrows |2008j Fernandez~| |2012[ i. The ex¬ 
plosion sets in shortly after M\ (Eq. 151 ) has fallen through 
the shock and well before the mass shell corresponding to 
M 4 + 0.3 has collapsed. High neutrino luminosities are 
maintained by high mass accretion rates and, after the onset 
of the explosion, by a contribution from the core emission 
(dashed lines in the luminosity panel of Fig. |2| that grows 
until roughly one second. The current models underestimate 
the surface luminosity of heavy-lepton neutrinos compared to 
more sophisticated simulations because of the chosen mod¬ 
est contraction of the inner boundary of the computational 
grid (which leads to underestimated temperatures in the accre¬ 
tion layer of the PNS) and because neutrino-pair production 
by nucleon-nucleon bremsstrahlung is not taken into account. 
We did not upgrade our treatment in this respect because v t , 
and v T are not of immediate relevance for our study since the 
explosion hinges exclusively on the heating by v e and v e . 

Replacing the inner core of the PNS by a contracting in¬ 
ner boundary of the computational mesh introduces a number 
of free parameters, whose settings allow one to achieve the 
desired accretion and neutrino-emission behavior as detailed 
above. On the one hand, our model contains parameters for 
the prescription of the contraction of the grid boundary, on the 
other hand there are parame ters for the simple high-density 
core model (cf. Ugliano et al.| [2012 and references therein). 
While the core-model parameters (T, f, R c (t), n) regulate the 
neutrino-emission evolution of the excised, high-density core 
of the PNS, the prescribed grid-boundary radius, R[b(t), gov¬ 
erns the settling of the hot accretion mantle of the PNS. Be¬ 
cause of partially compensating influences and dependences, 
not all of these parameters have a sensitive impact on the out¬ 
come of our study. Again, more relevant than a highly accu¬ 
rate description of individual components of the modeling is a 
reasonable reproduction of the overall properties of the accre¬ 
tion and neutrino emission history of the stalled SN shock and 
mass accumulating PNS. For example, the moderate increase 
of the mean neutrino energies with time and their regular hi¬ 
erarchy ((e Vf ) < (ev e ) < (%); Fig. El are not com patib le with 
the most sophisticated current models (Marek & Jank a|2009| 
Muller & Janka|20l4 see, e.g.,). They reflect our Choice ofa 


less extreme contraction of the 1.1 M 0 shell than found in sim¬ 
ulations with soft nuclear equations of state for the core mat¬ 
ter, where the PNS contracts more strongly and its accretion 
mantle heats up to higher temperatures at later post-bounce 
times (compare Scheck et al.||2006 and see the discussion by 
Pejcha & Thompson||2015|l. Our choice is motivated solely 


by numerical reasons (because of less stringent time-step con¬ 
straints), but it has no immediate drawbacks for our system¬ 
atic exploration of explosion conditions in large progenitor 
sets. Since neutrino-energy deposition depends on T v (e(r), the 
underestimated mean neutrino energies at late times can be 
compensated by higher neutrino luminosities L v of the PNS 
core. 

The neutrino emission from the PNS-core region is 
parametrized in accordance with basic physics constraints. 
This means that the total loss of electron-lepton number is 
compatible with the typical neutronization of the inner 1.1 M 0 
core, whose release of gravitational binding energy sati sfies 
energy conservation and virial theorem (see Ugliano et al. 


2012|. Correspondingly chosen boundary luminositiesThere- 
fore ensure a basically realistic deleptonization and cooling 
evolution of the PNS as a whole and of the accretion man¬ 
tle in particular, where much of the inner-boundary fluxes 
are absorbed and reprocessed. Again, a proper representa¬ 
tion of progenitor-dependent variations requires a reasonable 
description of the overall system behavior but does not need 
a very high sophistication of all individual components of the 
system. 

2.3.4. Calibration for low-ZAMS mass range 

Agreement with the constraints from SN 1987A employed 
in our work (i.e., the observed explosion energy, 56 Ni mass, 
total neutrino energy loss, and the duration of the neutrino sig¬ 
nal) can be achieved with different sets of values of the PNS 
core-model parameters. Using only one observed SN case the 
parameter set is underconstrained and the choice of suitable 
values is ambiguous. It is therefore not guaranteed that the 
calibration works equally well in the whole mass range of in¬ 
vestigated stellar models. 

In particular stars in the low ZAMS-mass regime (Mzams $ 
12-13 M 0 ) possess properties that are distinctly different from 
those of the adopted SN 1987A progenitors and in the mass 
neighborhood of these progenitors. Stars with Mzams ^ 12- 
13 M 0 are characterized by very small values of compactness 
(Eq. |T]», binding energy outside of the iron core and outside of 
M 4 (Eq. |2j, and mass derivative I-I 4 (Eq.[3|. The progenitor of 
SN 1054 giving birth to the Crab remnant is considered to be¬ 
long to this mass range, more specifically to have been a star 
with mass around 10 M Q ( |Nomoto et al.|[l982| |Smith||2013| 
Tominaga et al.|2013). Because of their structural similarities 
and distinctive differences compared to more massive stars, 
|Sukhbold et al.[ ( |2015[ ) call progenitors below roughly 12- 
13 M 0 “Crab-like” in contrast to stars above this mass limit, 
which they term “SN 1987-like”. 

Stars below ~10 M 0 were found to explode easily in self¬ 


consistent, sophisticated ID, 2D and 3D simulations 1 

Kitaura 

et al. 

2006; Janka et al. 

2008 

Fischer et al. 

2010 , 

Wanajo 

et al. 

|2011; Janka et al. 

2012 

Melson et a] 

. 2015a) with 


tie nickel production (< 0.01 M 0 ), in agreement with observa¬ 
tional properties concluded from detailed analyses of the Crab 
remnant (e.g.,|Yang & Chevalier|2015|). We therefore consider 
the results of these state-of-the-art SN models together with 
the empirical constraints for the Crab supernova as important 
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benchmark that should be reproduced by our approximate ID 
modeling of neutrino-powered explosions. 

The results of Ugliano et al. (|2012 1 revealed a problem 


in this respect, because they showed far more energetic ex 
plosions of stars in the low-mass domain than expected on 
grounds of the sophisticated simulations and from observa¬ 
tions of Crab. Obviously, the neutrino-source calibration used 
by Ugliano et al. ( |2012| > is not appropriate to reproduce “re¬ 


alistic” explosion conditions in stars with very dilute shells 
around the iron core. Instead, it leads to an overestimated 
power of the neutrino-driven wind and therefore overesti¬ 
mated explosion energies. In particular, the strong and ener¬ 
getic wind is in conflict with the short period of simultaneous 
postshock accretion and mass ejection after the onset of the 
explosion in <1OM 0 stars. Since the mass accretion rate is 
low and the duration of the accretion phase is limited by the 
fast shock expansion, the energetic importance of this phase 
is diminished by the small mass that is channelled through 


the neutrino-heating layer in convective flows (Kitaura et al. 

2006 Janka et al.J2008 

Wanajo et al. 2011 Janka et al.|]20l2f 

Melson et alTpOlSa). 

n order to account for these features 


found in the most refined simulations of low-mass stellar ex¬ 
plosions, the neutrino-driven wind power of our parametric 
models has to be reduced. 

We realize such a reduction of the wind power by decreas¬ 
ing the parameter which scales the compression work ex¬ 
erted on the inner (excised) core of the PNS by the overlying 
accretion mantle (cf. Eqs. 1-4 in Ugliano et al. [2012| , in pro¬ 
portionality to the compactness parameter ^ 1.75 b, which drops 
strongly for low-mass progenitors (see Eq. (3j). This proce¬ 
dure can be justified by the much lighter accretion layers of 
such stars, which implies less compression of the PNS core 
by the outer weight. Such a modification reduces the neutrino 
emission of the high-density core and therefore the mass out¬ 
flow in the early neutrino-driven wind. As a consequence, the 
explosion energy falls off towards the low-mass end of the in¬ 
vestigated progenitor sets. This can be seen in the upper panel 
of Fig. [3 which should be compared to the upper left panel of 
Fig. 5 in Ugliano et al. (2012). 

The £, scaling of Eq. Q is introduced as a quick fix in the 
course of this work and is a fairly ad hoc measure to cure the 
problem of overestimated explosion energies for low-mass SN 
progenitors. In|Sukhbold et al.|(j2015|l a different approach is 
taken, in which the final value of the core-radius parameter of 
the one-zone model describing the supernuclear PNS interior 
as neutrino source (see Ugliano et al. |2012| ) is modified. This 
procedure can directly be motivated by the contraction of the 
PNS found in self-consistent cooling simulations with micro¬ 
physical high-density equations of state and detailed neutrino 
transport. Mathematically, the modification of the core radius 
has a similar effect on the core-neutrino emission as the £ scal¬ 
ing employed here. While we refer the reader to |Sukhbold etj 
[ah ( f 201 5j ) for details, we emphasize that the consequences for 
the overall explosion behavior of the low-mass progenitors is 
very similar for both the £ reduction and the core-radius ad¬ 
justment applied by Sukhbold et al.|(|2015|). They lead to con¬ 
siderably lower explosion energies for Mzams £ 12 M 0 stars 
and a further drop of the explosion energies below —9.5 M 0 . 

As a drawback of this modification, the explosions of some 
of the low-mass progenitors between -10.5 M 0 and -12.5 M 0 
set in rather late (>1 s p.b., cf. Fig. [3|£] This, however, is ba¬ 
sically compatible with the tendency of relatively slow shock 


' In extremely rare cases one may even obtain failed explosions. 



Fig. 5. — Correlation of iron-core mass, Mp e , and compactness £ 5 for the 
investigated models of all progenitor series. Note that is taken from 
the pre-collapse model while f 1 5 is evaluated for a central density of 5 X 
10 10 gcrrT 3 . 


expansion and late explosions that are also found in sophisti¬ 
cated multi-dimensional simulations of such stars, which, in 
addition, reveal long-lasting phases of simultaneous accretion 
and mass ejection after the onset of the explosion (Muller 
et al. 2013 [Muller 2015| l. This extended accretion phase 
has only moderate consequences for the estimated remnant 
masses, because the mass accretion rate of these progenitors 
reaches a low level of <0.05-0.1 M 0 s _1 after a few 100ms 
post bounce, and some or even most of the accreted mass is 
re-ejected in the neutrino-driven wind. 

Besides providing information on the explosion energies 
and the onset times of the explosion (defined by the time 
the outgoing shock reaches 500 km), Fig. [3] also displays 
for exploding models the ejected masses o ~ (> Ni and iron- 
group tracer element, the baryonic and gravitational remnant 
masses, the fallback masses, and the total energies radiated 
by neutrinos. Overall, these results exhibit features very sim¬ 
ilar to those discussed in detail by Ugliano et al. ((2012]) for 
a different progenitor series and a different calibration model. 
We point out that the fallback masses in the low -mas s range 
of progenitors were overestimated by Ugliano et al. (12012 1 
due t o an error in the analysis (more discussion will follow in 
Sect. [375] an erratum on this aspect is in preparation.) 


3. RESULTS 

3.1. One and two-parameter classifications 

Figure [4] shows £ 2.5 versus ZAMS mass with BH forma¬ 
tion cases indicated by gray and explosions by red bars for 
the s2002 and s2014 series and all calibrations. The irregu¬ 
lar pattern found by Ugliano et al. (j2012|) for the s2002 pro¬ 
genitors is reproduced ancfappears similarly in the s2014 set. 
High compactness £ 2.5 exhibits a tendency to correlate with 
BHs. But also other parameters reflect this trend, for exam- 
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Fig. 6. — Correspondence of L v -M plane with critical neutrino luminosity L v , cr j t (M) (Ze/j) and x-y plane with separation curve y S ep(*) {right). In the left plot post¬ 
bounce evolution paths of successfully exploding models (white circles) and non-exploding models (black circles) are schematically indicated, corresponding to 
white and black circles for pre-collapse models in the right plot. Evolution paths of successful models cross the critical line at some point and the accretion ends 
after the explosion has taken off. In contrast, the tracks of failing cases never reach the critical conditions for launching the runaway expansion of the shock. 
The symbols in the left plot mark the “optimal point” relative to the critical curve that can be reached, corresponding to the stellar conditions described by the 
parameters (M 4 /i 4 ,/i 4 ) at the .9 = 4 location, which seems decisive for the success or failure of the explosion of a progenitor, because the accretion rate drops 
strongly outside. 


pie £ 15 , the iron-core mass M Fe (defined as the core where 

haus et al.| 2010 ; 

Hanke et al. 

20121 

Fernandez|2012j Pejcha 

E|i|A,> 46 )X/ > 0.5 for nuclei with mass numbers A,- and mass 

& Thompson 201 

2 Janka 201 

| 

ler & Janka 2015). M Fe 


fractions X ,), and the enclosed mass at the bottom of the O- 
buming shell. All three of them are tightly correlated, see 
Fig. Bias well as Fig. 4 of Ugliano et al. ( 2012| ). Also high val¬ 
ues of the binding energy Efim > M Fe ) outside of M Fe signals 
a tendency for BH formation, because this energy correlates 


with £>.5 (cf. Fig. 4 in Ugliano et af]2012). However, for none 
of these single parameters a sharp boundary value exists that 
discriminates between explosions and non-explosions. For all 
such choices of a parameter, the BH formation limit tends to 
vary (non-monotonically) with M Z a ms and in a broad interval 
of values either explosion or BH formation can happen. |Pe-| 
jcha & Thompson (|2015|) tried to optimize the choice of M 
for but even their best case achieved only 88% of correct 
predictions. Since in the cases of £>.5 and Efim > M Fe ), for 
example, the threshold value for BH formation tends to grow 
with higher ZAMS mass, one may hypothesize that a second 
parameter could improve the predictions. 

Placing the progenitors in a two-parameter space spanned 
by £ 1.5 and £ 2.5 or, equally good, M Fe and £ 2 . 5 , begins to show 
a cleaner separation of successful and failed explosions: SNe 
are obtained for small values of £ 2 . 5 , whereas BHs are formed 
for high values of £ 2 . 5 , but the value of this threshold increases 
with £15 and M Fe . For given f 1 .5 (or M Fe ) there is a value of 
£ 2.5 above which only BHs are formed. However, there is still 
a broad overlap region of mixed cases. 

This beginning separation can be understood in view of 
the theoretical background of the neutrino-driven mechanism, 
where the expansion of the SN shock is obstructed by the 
ram pressure of infalling stellar-core matter and shock ex¬ 
pansion is pushed by neutrino-energy deposition behind the 
shock. For neutrino luminosities above a critical threshold 
Ty.crit (M), which depends on the mass-accretion rate M of the 
shock, shock runaway and explosion are triggered by neu 
trino heating (see Fig. P" 

IJankal 20C 


left panel, and, e.g.. 


Burrows & 


, j Murphy & Burrows||2008f|Nord- 


_j of the 

PNS as accretor, which determines the strength of the grav¬ 
itational potential and the size of the neutrino luminosities. 
Such a dependence can be concluded from the proportional¬ 
ity L v oc RyTy, where R v is the largely progenitor-independent 
neutrinosphere radius and the neutrinospheric temperature 
T v inc reases roughly linearly with M ns (see Muller & Janka] 
|2014| >. On the other hand, the long-time mass-accretion rate 
of the PNS grows with £ 2 . 5 . which is higher for denser stel¬ 
lar cores. For each PNS mass explosions become impossible 
above a certain value of M or £ 2.5 • 

3.2. Two-parameter classification based on the theoretical 
concept of the neutrino-driven mechanism 

If the initial mass cut at the onset of the explosion devel¬ 
ops at an enclosed mass M = m(r) of the progenitor, we can 
choose M = m(r) as a suitable proxy of the initial PNS mass, 
M ns . A rough measure of the mass-accretion rate M by the 
stalled shock around the onset of the explosion is then given 
by the mass-gradient m'(r) = dm(r)/dr = A::r 2 p(r) at the cor¬ 
responding radius r. This is the case because m'(r) can be 
directly linked to the free-fall accretion rate of matter collaps¬ 
ing into the shock from initial radius r according to 


^ dm dm(r) / dfff(r) \ 
dfff dr \ dr / 
2 m' (r) 


-1 


hf[(3/r) - m'(r)/m(r )] 


2 r 

- m (r), 

3 f ff 


(5) 



core. 


Following the critical-luminosity concept now points the 
way to further improvements towards a classification scheme 

































































12 


Ertl et al. 


of explosion conditions: The L v -M dependence of the 
neutrino-driven mechanism suggests that the explodability of 
the progenitors may be classified by the parameters M = m(r) 
and M oc m'{r), because the accretion luminosity L‘ [ , cc oc 
GM ns M/R ns oc Mm'(r) accounts for a major fraction of the 
neutrino luminosity of the PNS at the time of shock revival 
( [Miiller & Janka||2014| Muller & Janka||20 1 5| >, and, in par¬ 
ticular, it is the part of the neutrino emission that reflects the 
main progenitor dependence. It is important to note that the 
time-evolving NS and neutrinospheric radii, R ns ~ R v , are 
nearly the same for different progenitors and only weakly time 
dependent when the explosions take place rather late after 
bounce. This is true for our simulations (where t exp > 0.5 s 
with few exceptions; Fig. [3]) as well a s for self-consistent, so 
phisticated models (see, e.g.. Fig. 3 in Muller & Janka|20l4 


In both cases the spread of the NS radii and their evolution be¬ 
tween ~0.4 s and 1 s after bounce accounts for less than 25% 
variation aro und a n average value of all investigated mod¬ 
els. In Sect. |3.6| we will come back to this argument and 
give reasons why the NS radius has little influence on the re¬ 
sults discussed in this work. Moreover, the neutrino loss from 
the low-entropy, degenerate PNS core, whose properties are 
determined by the incompressibility of supranuclear matter, 
should exhibit a progenitor dependence mostly through the 
different weight of the surrounding accretion mantle, whose 
growth depends on M. Such a connection is expressed by the 
terms depending on ;n acc and m acc in the neutrino luminos¬ 


ity of the high-density PNS core in Eq. (4) of Ugliano et al. 
(2012l. We therefore hypothesize, and demonstrate below, a 


correspondence of the L v -M space and the Mm' -in' parameter 
plane and expect that the critical luminosity curve Ly C „ t (M) 
maps to a curve separating BH formation and successful ex¬ 
plosions in the Mm'-m' plan^] 

In our simulations neutrino-driven explosions set in around 
the time or shortly after the moment when infalling matter ar¬ 
riving at the shock possesses an entropy s ~ 4. We therefore 
choose M 4 (Eq.|2| as our proxy of the PNS mass, M 4 oc M ns , 
andyi /4 = m'(r)[M 0 /(lOOOkm )] _1 | s= 4 (Eq. as corresponding 
measure of the mass-accretion rate at this time, p 4 oc M. The 
product Mm' is therefore repesented by M 4 J u 4 . Tests showed 
that replacing M 4 by the iron-core mass, M Fe , is similarly 
good and yields results of nearly the same quality in the analy¬ 
sis following below (which points to an underlying correlation 
between M 4 and M Fe ). In practice, we evaluate Eq. 0 for /i 4 
by the average mass-gradient of the progenitor just outside of 
5 = 4 according to 


P4 = 


A m/M 0 


Ar/lOOOkm 


s =4 


(M 4 + A m/M 0 ) - M 4 


[r(M 4 + A m/M 0 ) - r(s = 4)] /1000 km 


( 6 ) 


with Am = 0.3 M G yielding optimal results according to tests 
with varied mass intervals A in. With the parameters M 4 and fj 4 
picked, our imagined mapping between critical conditions in 
the Ly-M and Mm' -in' spaces transforms into such a mapping 

6 Since the shock revival is determined by neutrino heating, which depends 
on Ly(£y), and since the average sq uared neutrino energy ( e} : ) oc 7V; oc /Vf/, 
increases with time and PNS mass I Muller & Janka 2014 1 , Muller & Janka 
(2015) discuss the critical condition for shock revival in terms of L v { 4 ) as 
a function of M ns M. This suggests that an alternative choice of parameters 
could be M 3 m' and Mm' instead of Mm' and m', respectively. Our results 
demonstrate that the basic physics is already captured by the Mm' -m' depen¬ 
dence. 


relation between the L y -M and M 4 // 4-//4 planes as illustrated 

by Fig. [6] 

Figure [7] demonstrates the strong correlation of the mass 
accretion rate M - dm/dt with the parameter p 4 as given 
by Eq. ([ 6 ]) (panel c) as well as the tight correlations between 
/V/ 4/./4 and the sum of v e and v e luminosities (panel a) and the 
summed product of the luminosities and mean squared ener¬ 
gies of v e and v e (panel b). It is important to note that the non- 
stationarity of the conditions requires us to average the quanti¬ 
ties plotted on the abscissas over time from the moment when 
the 5 = 4 interface passes through the shock until the models 
either explode (i.e., the shock radius expands beyond 500 km; 
open circles) or until the mass shell (M 4 + 0.3) M Q has fallen 
through the shock, which sets an endpoint to the time interval 
within which explosions are obtained (non-exploding cases 
marked by filled circles). The time averaging is needed not 
only because of evolutionary changes of the preshock mass- 
accretion rate (as determined by the progenitor structure) and 
corresponding evolutionary trends of the emitted neutrino lu¬ 
minosities and mean energies. The averaging is necessary, in 
particular, because the majority of our models develops large- 
amplitude shock oscillations after the accretion of the 5 = 4 
interface, which leads to quasi-periodic variations of the neu¬ 
trino emission properties with more or less pronounced, grow¬ 
ing amplitudes (see the examples in Fig. |2|. Panel d demon¬ 
strates that exploding models (open circles) exceed a value of 
unity for the ratio of advection time scale, t a dv> to heating time 
scale, fheat, in the gain layer, which was considered as a use- 


vious works (e.g., Janka & Keil 1998; Thompson 2000 

Janka 

et al. 2001; Thompson et al. 2005, 

Buras et al.]2006; 

Marek 

& Janka 2009; 

Miiller et al. 2012] 

d\ IMiiller & Janka [201 5[ 


gion towards low mass-accretion rates (as visible in pane ls a- 
c, too), which confirms our observation reported in Sect. 3.1 
In contrast, non-exploding models cluster, clearly separated, 
in the left, upper area of panel d, where 4 dv/ 4 eat £ 1 and the 
mass-accretion rate tends to be higher. For the calculation of 
the t ime scales we follow the definitions previously used by, 
e.g.,|Buras et al.|(|2006|l,|Marek & Janka|(|2009), Muller et al.| 
d2012bK , |Muller & Janka|fl2015| ): 


r R* 


fiieat — 


I" (e + 0)pdV 

jR e 

_ r Rs 1 , 

^adv — I , . dr. 

JR „ I v r 


\ r>R & \ _1 

I q v p dv 
Jr „ 


(7) 

( 8 ) 


Here, the volume and radius integrals are performed over the 
gain layer between gain radius R„ and shock radius R s . e is 
the sum of the specific kinetic and internal energies, <t> the 
(Newtonian) gravitational potential, p the density, q Y the net 
heating rate per unit of mass, and v r the velocity of the flow. 
Again, because of the variations of the diagnostic quantities 
associated with the time evolution of the collapsing star and 
the oscillations of the gain layer, the mass-accretion rate and 
time-scale ratio are time-averaged between the moment when 
the 5 = 4 interface passes the shock until either 300 ms later 
or until the model explodes (shock radius exceeding 500 kmn 
Panels e (and f) lend support to the concept of a critical thresm 
old luminosity in the L v -M (and L v {el)-M) space mentioned 

7 We tested intervals ranging from 100 ms to 600 ms and observed the same 
trends for all choices. 
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Fig. 7.— End points of the postbounce evolution of exploding (open circles) and non-exploding models (filled circles) in planes spanned by various pairs of 
parameters. The data correspond to results of the s2002 model series with wl8.0 calibration. All symbols represent time-averaged locations because of strong 
temporal variations of the postshock accretion layer. Panels a-c display correlations of our dimensionless progenitor parameters p 4 and M 4//4 with time-averaged 
values of the preshock mass accretion rate, M = dm/dt, and the characteristic neutrino-emission properties, respectively, as obtained in our simulations and 
measured at 500 km. The abscissas of panels a and b show the summed luminosities of v e and v e , L Ve + Ly e , and the summed products of luminosities and mean 
squared energies, Ly e (€* e )+Ly e (e£) f of both neutrino species, respectively. For the exploding models the time averaging is performed from the arrival of the .9 = 4 
interface at the shock until the explosion sets in (defined by the shock radius reaching 500 km), whereas the averages for non-exploding models cover the time 
from the 5 = 4 interface passing the shock until 0.3 M Q of overlying material have been accreted by the shock. Panel d displays the separation of exploding and 
non-exploding models in the plane spanned by the mass-accretion rate and the ratio of advection to heating time scale. Panels e andf demonstrate this separation 
in the planes spanned by M and L Ve + Ly e or L Ve (€y e ) + Ly e (e? ), respectively. The time averages of the quantities in panels d-f are computed from the passage 
of the 5 = 4 interface through the shock until 300 ms later for non-exploding models or until the onset of the explosion otherwise. Gray shading in panels d-f 
indicates the regions where explosions fail. (No exact boundary curves are determined for the cases of panels e andf.) 
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above. The two plots show a separation of exploding (open 
circles) and non-exploding (filled circles) models in a plane 
spanned by the time-averaged values of the preshock mass ac¬ 
cretion rate on the one hand and the sum of v e and v e luminosi¬ 
ties (panel e) or the summed product of v e and v e luminosities 
times their mean squared energies on the other hand (panel f). 
Again the same time-averaging procedure as for panel d is 
applied. After the s = 4 progenitor shell passes the shock, 
the time-averaged conditions of the exploding models reach 
the lower halfs of these panels, whereas the time-averaged 
properties of the non-exploding models define the positions 
of these unsuccessful explosions in the upper halfs. A sep¬ 
aration appears that can be imagined to resemble the critical 
luminosity curve L v , cr j t (M) sketched in the left panel of Fig.pl 

Because of the strong time dependence of the postshock 
conditions and of the neutrino emission during the phase of 
dynamical shock expansion and contraction after the accre¬ 
tion of the 5 = 4 interface, it is very difficult to exactly de¬ 
termine the critical luminosity curve that captures the physics 
of our exploding ID simulations. Although we do not con¬ 
sider such an effort as hopeless if the governing parameters are 
carefully taken into account (cf., e.g., the disc ussions in |Pe-| 
jcha & Thompsonj 2015 ||Miiller & Janka| 2015 ]l and their time 
variations are suitably averaged, we think that the clear sepa¬ 
ration of successful and failed explosions visible in panels e 
and f of Fig. IT] provides proper and sufficient support for the 
notion of such a critical curve (or, more general: condition) 
in the L VCI i t -M space. Figure [6] illustrates our imagined rela¬ 
tion between the evolution tracks of collapsing stellar cores 
that explode or fail to explode in this L v -M space on the left 
side and the locations of SN-producing and BH-forming pro¬ 
genitors in the M 4/24-/14 plane on the right side. The sketched 
evolution paths are guided by our results in panels e and f of 
Fig. [ 7 ] In the following section we will demonstrate that ex¬ 
ploding and non-exploding simulations indeed separate in the 
M4/14-/14 parameter space. 

3.3. Separation line of exploding and non-exploding 
progenitors 

The existence of a separation line between BH forming and 
SN producing progenitors in the M4/14-/24 plane is demon¬ 
strated by Fig. [8] which shows the positions of the progenitors 
for all investigated model series in this two-dimensional pa¬ 
rameter space. For all five calibrations successful explosions 
are marked by colored symbols, whereas BH formation is in¬ 
dicated by black symbols. The regions of failed explosions are 
underlaid by gray. They are bounded by straight lines with fit 
functions as indicated in the panels of Fig. [8] 

JsepM = k\ ■ X + h , (9) 

where x = A7 4 // 4 and y = p4 are dimensionless variables with 
M4 in solar masses and // 4 computed by Eq. j6). The values 
of the dimensionless coefficients k\ and k 2 as listed in Table [ 2 ] 
are determined by minimizing the numbers of outliers. 

The stellar models of all progenitor sets populate a narrow 
strip in the x-y plane of Fig. IS] left panels. BH formation 
cases are located in the upper left part of the x-y plane. The 
inclination of the separation line implies that the explosion 
limit in terms of P4 depends on the value of the M 4 /i 4 and 
therefore a single parameter would fail to predict the right be¬ 
havior in a large number of cases. Denser cores outside of M4 
with high mass-accretion rates (larger P4) prevent explosions 
above some limiting value. This limit grows for more massive 


TABLE 2 

Bh-sn separation curves for all calibration models 


Calibration Model 

*l a 

*2 a 

M 4 b 



sl9.8 (2002) 

0.274 

0.0470 

1.529 

0.0662 

0.101 

wl5.0 c 

0.225 

0.0495 

1.318 

0.0176 

0.023 

wl8.0 

0.283 

0.0430 

1.472 

0.0530 

0.078 

w20.0 

0.284 

0.0393 

1.616 

0.0469 

0.076 

n20.0 

0.194 

0.0580 

1.679 

0.0441 

0.074 


a Fit parameters of separation curve (Kq. [7] when x and y are measured for a 
central stellar density of 5 X 10 10 g cm~ 3 . 
b Measured for a central stellar density of 5 X 10 lo gcm~ 3 . 
c M\ and p 4 measured roughly at core bounce, because pre-collapse data are 
not available. 

cores and thus higher M4 because larger mass-accretion rates 
do not only hamper shock expansion by higher ram pressure 
but larger core masses and bigger accretion rates also corre¬ 
late with an increase of the neutrino luminosity of the PNS 
as expressed by our parameter M4P4. The evolution tracks of 
successful explosion cases in the left panel of Fig. [6] indicate 
that for higher L v oc M ns M the explosion threshold, L vcnl (M), 
can be reached for larger values of M. 

Explosions are supported by the combination of a massive 
PNS, which is associated with a high neutrino luminosity 
from the cooling of the accretion mantle, on the one hand, 
and a rapid decline of the accretion rate, which leads to de¬ 
creasing ram pressure, on the other hand. A high value of M4 
combined with a low value of P4 is therefore favorable for an 
explosion because a high accretion luminosity (due to a high 
accretor mass M ns ~ M 4 ) comes together with a low mass 
accretion rate (and thus low ram pressure and low binding 
energy) exterior to the 5 = 4 interface (cf. Fig[9]». Such condi¬ 
tions are met, and explosions occur readily, when the entropy 
step at the 5 = 4 location is big, because a high entropy value 
outside of M4 correlates with low densities and a low accre¬ 
tion rate. M4 is usually the base of the oxygen shell and a 
place where the entropy changes discontinuously causing (or 
resulting from) a sudden decrease in density due to burning 
there. This translates into an abrupt decrease in M when the 
mass M 4 accretes. Figure 14 of |Sukhbold & Woosley| ( j2014) i 
shows a strong correlation between compactness £2.5 and lo¬ 
cation of the oxygen shell. The decrease of the mass accretion 
rate is abrupt only if the entropy change is steep with mass, 
for which // 4 at M 4 is a relevant measure. 

Progenitors with Mzams 5! 22 M G that are harder to explode 
often have relatively small values of M4 and an entropy ledge 
above 5 = 4 on a lower level than the entropy reached in 
more easily exploding stars. The lower neutrino luminosity 
associated with the smaller accretor mass in combination with 
the higher ram pressure can prohibit shock expansion in many 
of these cases. Corresponding to the relatively small values of 
M4 and relatively higher densities outside of this mass, these 
cases stick out from their neighboring stars with respect to the 
binding energy of overlying material, namely, non-exploding 
models in almost all cases are characterized by local maxima 
of E hA = Eb(m/M Q > M4) (see Fig. [9]. 

In view of this insight it is not astonishing that exploding 
and non-exploding progenitors can be seen to start separat¬ 
ing from each other in the two-parameter space spanned by 
M 4 and the average entropy value (5)4 just outside of M4 
(Fig. [TO] !. Averaging 5 over the mass interval [M 4 , M 4 + 0.5] 
turns out to yield the best results. Exploding models clus¬ 
ter towards the side of high (5)4 and low M4, while failures 
are found preferentially for low values of (5)4. The threshold 
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Fig. 8. — Separation curves between BH formation (gray region, black symbols) and SN explosions (white region, colored symbols) for all calibrations in the 
plane of parameters x = V/4//4 and y = /14 (zooms in right panels). Note that the left panels do not show roughly two dozen BH-forming models of the u2002 
series, which populate the x-range between 0.5 and 0.62 and are off the displayed scale. Different symbols and colors correspond to the different progenitor sets. 
The locations of the calibration models are also indicated in the left panels by crossing blue lines. 


for success tends to grow with M 4 . However, there is still a 
broad band where both types of outcomes overlap. The disen¬ 
tanglement of SNe and BH-formation events is clearly better 
achieved by the parameter set of /V/ 4 /j 4 and // 4 , which, in addi¬ 
tion, applies correctly not only for stars with Mzams > 15 M 0 
but also for progenitors with lower masses. 

3.4. Stellar outliers 

Out of 621 simulated stellar models for the sl9.8, wl5.0, 
wl8.0, w20.0, and n20.0 calibrations only 9, 14, 16, 11, and 9 
models, respectively, do not follow the behavior predicted by 
their locations on the one or the other side of the separation 
line in the M 4 / 24-/24 plane (see the zooms in the right column 
of Fig. [ 8 ]). But most of these cases lie very close to the bound¬ 
ary curve and their explosion or non-explosion can be affected 
by fine details and will certainly depend on multi-dimensional 
effects. A small sample of outliers is farther away from the 
boundary line. The w20.0 calibration is the weakest driver of 
neutrino-powered explosions in our set and tends to yield the 
largest number of such more extreme outliers. 


These cases possess unusual structural features that influ¬ 
ence their readiness to explode. On the non-exploding side 
of the separation line, model s20.8 of the s2014 series with 
(M 4 /i 4 ,// 4 ) ~ (0.142,0.0981) is one example of a progenitor 
that blows up with all calibrations except w 20 . 0 , although it 
is predicted to fail (see Fig. [ 8 ]). In contrast, its mass-neighbor 
s20.9 with (M 4 / 4 ,/ 2 4 ) * (0T23,0.085) as well as its close 
neighbor in the M 4 / 24-/24 space, sl5.8 of the s2002 series with 
(M 4 // 4 ,/ 2 4 ) ~ (0.140,0.096), both form BHs as expected. The 
structure of these pre-supernova models in the s — 4 region 
is very similar with M 4 = 1.45,1.45,1.46 for s20.8, s20.9, 
sl5.8, respectively. Although s20.8 reaches a lower entropy 
level outside of s — 4 than the other two cases and therefore 
is also predicted to fail, its explosion becomes possible when 
the next entropy step at an enclosed mass of 1.77 M Q reaches 
the shock. This step is slightly farther out (at 1.78 M 0 ) in the 
s20.9 case and comes much later (at ~1.9M 0 ) in the s 15.8 
model. Both the earlier entropy jump and the lower preced¬ 
ing entropy level enable the explosion of s 20 . 8 , because the 
associated higher density maintains a higher mass-accretion 
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Fig. 9. — Iron-core masses Mp e {top panel), exterior binding energies (£b,Fe = Eb{m > Mp e ); second panel), normalized masses M 4 {third panel), exterior 
binding energies (£b ,4 = E^ui/Mq > M 4 ); fourth panel), and (fifth panel), for models of the s2014 series, the supplementary low-mass progenitors with 
Mzams < 15 M g , and models with Mzams > 30 M© from the s2007 series. A black vertical line marks the boundary between the two progenitor sets. Red 
bars indicate exploding cases and gray bars non-exploding ones. All quantities are measured when the central density of the collapsing stellar iron core is 
5 x 10 10 gcm -3 . The upper five panels correspond to the wl8.0 calibration, the lower five panels to the n20.0 calibration. Mass and parameter values of 
the calibration models are indicated by vertical and horizontal blue lines, respectively. In the region of Mzams ^ 22 M© non-exploding cases, with very few 
exceptions, correlate with local minima of M 4 and pronounced local maxima of p 4 and £^, 4 - A high value of M 4 combined with a low value of is typically 
supportive for an explosion because a high accretion luminosity (due to a high accretor mass M ns « M 4 ) comes together with a low mass accretion rate (and thus 
low ram pressure and low binding energy) exterior to the 5 = 4 interface. The iron-core masses and their exterior binding energies show a similar tendency, but 
significantly less pronounced. 
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Fig. 10.— Two-parameter plane spanned by M 4 = M(s = 4) and mean 
entropy of overlying matter, < 5 ) 4 , averaged over a mass intervall of 0.5 M©. 
The locations of progenitors from the s2014 series and the supplementary 
low-mass set with Mzams < 15 M© are marked by crosses. Black crosses 
indicate BH formation cases and tend to concentrate towards the lower and 
right half of the panel, i.e., towards low values of < 5)4 for given M 4 . Orange 
crosses mean successful SN explosions of the s2014 models with the wl8.0 
calibration, and red crosses are explosions for the Mzams <15 M© progen¬ 
itors. While high entropies outside of the s = 4 location signal a tendency 
of success for the stars of the s2014 set (although the separation from BH- 
formation cases is not sharp), the successfully exploding Mzams < 15M© 
models mix completely with BH forming events. 


rate and therefore higher neutrino luminosity until the entropy 
jump at 1.77 M 0 falls into the shock. The abnormal structure 
of the progenitor therefore prevents that the explosion behav¬ 
ior is correctly captured by our two-parameter criterion for the 
explodability. 

On the exploding side, model sl5.3 of the s2014 series with 
pf) * (0.146,0.0797) is expected to blow up accord¬ 
ing to the two-parameter criterion, but does not do so for all 
calibrations (Fig. [8]). Similarly, sl5.0 of the s2007 series with 
(M^pn,pf) » (0.137,0.0749) fails with the wl5.0 and w20.0 
calibrations although success is predicted. We compare their 
structure with the nearby successful cases of s25.4 (s2014 se¬ 
ries, {M^p^,pf) =3 (0.150,0.0820)), and s25.2, s25.5 (both 
from the s2014 series), and s25.8 (s2002 series), all of which 
group around as (0.143,0.0783). The success¬ 

fully exploding models all have similar entropy and density 
structures, namely fairly low entropies (s < 3) and therefore 
high densities up to 1.81-1.82 M 0 , where the entropy jumps 
to s > 6. The high mass-accretion rate leads to an early arrival 
of the s = 4 interface at the shock (~300 ms after bounce) and 
high accretion luminosity. Together with the strong decline 
of the accretion rate afterwards this fosters the explosion. In 
contrast, the two models that blow up less easily have higher 
entropies and lower densities so that the s = 4 mass shells (at 
~1.8 M 0 in sl5.0 and at ~ 1.82 M 0 in the s 15.3) arrive at the 
shock much later (at ~680ms and ~830ms post bounce, re¬ 
spectively), at which time accretion contributes less neutrino 
luminosity. Moreover, both models have a pronounced en¬ 
tropy ledge with a width of ~O.O5M 0 (sl5.0) and ~0.08 M 0 
(sl5.3) before the entropy rises above s ~ 5. This ledge is 
much narrower than in the majority of non-exploding mod¬ 
els, where it stretches across typically 0.3 M 0 or more. The 
continued, relatively high accretion rate prohibits shock ex¬ 
pansion and explosion. This is obvious from the fact that 
model sl5.0 with the less extended entropy ledge exhibits a 



X = M 4 yU.., 


Fig. 11.— ZAMS masses, iron-core masses, Mp e , and binding energies 
outside of the iron core, E\^(m > Mp e ), (from top to bottom ) of the s2014 
series and the supplementary low-mass progenitors with Mzams < 15M© 
in the x-y parameter plane. The solid and dashed lines mark the separation 
curves y S ep(*) for the wl 8.0 and n 20.0 calibration models (different symbols 
as given in the legend), respectively. 

stronger tendency to explode and for some calibrations indeed 
does, whereas sl5.3 with the wider ledge fails for all cali¬ 
brations. Our diagnostic parameter p\ to measure the mass 
derivative in an interval of Am = 0.3 M 0 , however, is domi¬ 
nated by the high-entropy level (low-density region) above the 
ledge and therefore underestimates the mass-accretion rate in 
the ledge domain, which is relevant for describing the explo¬ 
sion conditions. Again the abnormal structure of the si5.0 
and s 15.3 progenitors prevents our two-parameter classifica¬ 
tion from correctly describing the explosion behavior of these 
models. 

3.5. Systematics of progenitor and explosion properties in the 
two-parameter plane 

In Fig. [TT| colored symbols show the positions of the pro¬ 
genitors of the s2014 series and those of the supplementary 
low-mass models with M z ams <15 M 0 in the x-y-plane rel¬ 
ative to the separation lines y scp (x) of exploding and non¬ 
exploding cases. In the upper panel the color coding corre¬ 
sponds to Mzams> in the middle panel to the iron-core mass, 
Mp e , and in the bottom panel to the binding energy of matter 
outside of the iron core. Mp e is taken to be the value pro¬ 
vided by the stellar progenitor model at the start of the col¬ 
lapse simulation in order to avoid misestimation associated 
with our simplified nuclear burning network and with inaccu¬ 
racies from the initial mapping of the progenitor data. Since 
we use the pressure profile of the progenitor model instead 
of the temperature profile, slight differences of the derived 
temperatures can affect the temperature-sensitive shell burn¬ 
ing and thus the growth of the iron-core mass. 

While low-mass progenitors with small iron cores and low 
binding energies populate the region towards the lower left 
corner with significant distance to the separation curve, stars 
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Fig. 12.— Explosion energies (1 B = 1 bethe = 10 51 erg), post-bounce explosion times, gravitational neutron-star masses (M ns g = A7, ls y - E vtot /c 2 ), ejected 
iron-group material (i.e., 56 Ni plus tracer masses), and fallback masses (from top to bottom) of the s2014 series and the supplementary low-mass progenitors 
with Mzams < 15 M 0 for calibration models wl8.0 (left) and n20.0 (right) in the x-y parameter plane. Black crosses correspond to BH formation cases, colored 
crosses to successful explosions. In the middle and bottom panels, the blue (partly overlapping) symbols correspond to fallback SNe with estimated BH masses 
(baryonic masses in parentheses) and fallback masses as listed in the legends. The horizonal and vertical lines mark the locations of the calibration models with 
the colors corresponding to the values of the displayed quantities. 


above 20 M 0 with bigger iron cores and high binding energies 
can be mostly found well above the separation curve. How¬ 
ever, there are quite a number of intermediate-mass progeni¬ 
tors above the line and higher-mass cases below. In particular, 
a lot of stars with masses between -25 M 0 and 30 M 0 cluster 
around y sep (x) in the x — 0.13—0.15 region. These stars are 
characterized by Mp e ~ 1.4—1.5 M 0 and high exterior binding 
energies. Some of them explode but most fail (cf. Fig. [9|. 
The ones that group on the unsuccessful side are mostly cases 
with smaller iron cores, whose neutrino luminosity is insuffi¬ 
cient to create enough power of neutrino heating to overcome 
the ram pressure of the massive infall. 

Figure fl2| displays the BH-formation cases of the s2014 
series without associated SNe by black crosses in the x-y- 
plane. Successful SN explosions of this series plus additional 
Mzams < 15 M q progenitors are shown by color-coded sym¬ 
bols, which represent, from top to bottom, the final explosion 


energy (£e Xp , with the binding energy of the whole progen¬ 
itor taken into account), the explosion time (f exp , measured 
by the time the outgoing shock reaches 500 km), the gravita¬ 
tional mass of the remnant (with fallback taken into account), 
the ejected mass of 56 Ni plus tracer element (see Sect.fZj fall¬ 
back also taken into account), and the fallback mass. The left 
plot shows the results of our wl8.0 calibration, the right plot 
for the n20.0 calibration. The gravitational mass of the NS 
remnant, M ns g , is estimated from the baryonic mass, M nSj b, 
by subtracting the rest-mass equivalent of the total neutrino 
energy carried away in our simulations: 


ns.g — 'H ns h 


1 

— Fy tot • 


( 10 ) 


Our estimates of NS binding energies, E ns _\, — E Vttot , are 

\\\ 1 989|> fit of£j is Y h = 
1.5 x lO 53 (M ns , g /M 0 ) 2 erg = 0.084 M 0 c 2 (M ns>g /M 0 ) 2 . Blue 


roughly compatible with 


Fattimer & Yahil 
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symbols in the middle and bottom panels mark fallback BH 
formation cases, for which gravitational and baryonic masses 
(the latter in parentheses) are listed in the panels. We ignore 
neutrino-energy losses during fallback accretion for both NS 
and BH-forming remnants. 

Progenitors in Fig. fl2| that lie very close to the separation 
line tend to produce weaker explosions that set in later than 
those of progenitors with a somewhat greater distance from 
the line. Moreover, there is a tendency of more massive NSs 
to be produced higher up along the x-y-band where the pro¬ 
genitors cluster, i.e., bigger NS masses are made at higher 
values of x = M 4 / 44 . Also the largest ejecta masses of 56 Ni 
and 56 Ni plus tracer are found towards the right side of the 
displayed progenitor band just below the boundary of the BH 
formation region. 

Fallback masses tend to decrease towards the lower left 
corner of the x-v-plane in Fig. [T2] far away from the sepa¬ 
ration curve, where predominantly low-mass progenitors are 
located, besides five progenitors around 20 M Q , which lie 
in this region because they have extremely low values of 
E\,(m > M 4 ), see Fig. [9] and exceptionally small values of 
/r 4 (Fig. Hi, and develop fast and strong explosions with 
small NS masses, large masses of ejected 56 Ni plus tracer, and 
very little fallback. Closer to the separation line the fallback 
masses are higher, but for successfully exploding models they 
exceed ~0.05 M Q only in a few special cases of fallback SNe 
(cf. Fig. [3]t. where the fallback mass can amount up to several 
solar masses. We point out that the fallback masses in par¬ 
ticular of stars below ~20 M Q were massively overestimated 
by Ugliano et ah (|2012j). The reason was an erroneous in¬ 
terpretation of the outward reflection of reverse-shock accel¬ 
erated matter as a numerical artifact connected to the use of 
the condition at the inner grid boundary. The reverse shock, 
which forms when the SN shock passes the He/H interface, 
travels backward through the ejecta and decelerates the out¬ 
ward moving matter to initially negative velocities. This in¬ 
ward flow of stellar material, however, is slowed down and 
reflected back outward by the large negative pressure gradient 
that builds up in the reverse-shock heated inner region. With 
this outward reflection, which is a true physical phenomenon 
and not a boundary artifact, the matter that ultimately can be 
accreted by the NS is diminished to typically between some 
10~ 4 M 0 and some 10 2 M Q of early fallback (see Fig. [ 3 }. 

In a handful of high-mass s2014 progenitors (s27~2 and 
s27.3 for the wl8.0 calibation and s27.4, s29.0, s29.1, s29.2, 
and s29.6 for the n20.0 set) the SN explosion is unable to 
unbind a large fraction of the star so that fallback of more 
than a solar mass of stellar matter is likely to push the NS 
beyond the BH formation limit. Such fallback SN cases clus¬ 
ter in the vicinity of x ~ 0.13-0.14 and y ~ 0.080-0.081 on 
the explosion side of the separation line. In the ZAMS-mass 
sequence they lie at interfaces between mass intervals of suc¬ 
cessfully exploding and non-exploding models, or they appear 
isolated in BH-formation regions of the ZAMS-mass space 
(see Figs. Bland|9). Their fallback masses and estimated BH 
masses arelistedin the corresponding panels of Fig. 12 Nat¬ 
urally, they stick out also by their extremely low ejecta masses 
of 56 Ni and tracer elements, late explosion times (around one 
second post bounce or later), and relatively low explosion en¬ 
ergies (~0.3-0.5 B). 


3.6. Influence of the neutron star radius 

Basically, the accretion luminosity, which is given by Lf c oc 
GM ns M/R ns , does not only depend on the PNS mass, M ns , 
and the mass-accretion rate, M, but also on the PNS radius, 
R ns . One may wonder whether our two-parameter criterion 
is able to capture the essential physics although we disregard 
the radius dependence when using x = M 4 /r 4 as a proxy of the 
accretion luminosity. 

For this reason, we tested a redefinition of x = M 4/14 by 
including a factor (R n scorns’ i- e ” usin 8 * s 
insteaqj Doing so, we found essentially no relevant effects 
on the location of the boundary curve. In fact, the sepa¬ 
ration of exploding and non-exploding models in the x-y- 
plane is even slightly improved compared to the x-v-plane, 
because (/?ns)J(io ms f° r most non-exploding models is larger 
than for the far majority of exploding ones. As a consequence, 
the non-exploding cases tend to be shifted to the left away 
from the boundary line, whereas most of the exploding cases 
are shifted to the right, also increasing their distances to the 
boundary line. This trend leads to a marginally clearer dis¬ 
entanglement of both model groups near the border between 
explosion and non-explosion regions. A subset of the (any¬ 
way few and marginal) outliers can thus move to the correct 
side, while very few cases can become new, marginal outliers. 
It might therefore even be possible to improve the success 
rate for the classification of explodability by a corresponding 
(minor) relocation of the boundary curve. The improvement, 
however, is not significant enough to justify the introduction 
of an additional parameter into our two-parameter criterion 
in the form of (R ns ) 3 ooms> which has the disadvantage of not 
being based on progenitor properties and whose exact, case- 
dependent value cannot be predicted by simple arguments. 

Because the line separating exploding and non-exploding 
models did not change in our test, the criterion advocated in 
this paper captures the basic physics, within the limitations 
of the modeling. While we report here this marginal sensitiv¬ 
ity of our two-parameter criterion to the NS radius as a result 
of the present study, a finally conclusive assessment of this 
question would require to repeat our set of model calculations 
for different prescriptions of the time-dependent contraction 
of the inner boundary of our computational grid. The cho¬ 
sen functional behavior of this boundary movement with time 
determines how the proto-neutron star contraction proceeds 
during the crucial phase of shock revival. In order to avoid 
overly severe numerical time-step constraints, which can be¬ 
come a serious handicap for our long-time simulations with 
explicit neutrino tr anspor t over typically 20 seconds, we fol¬ 
low Ugliano et al.] P012jl in using a relatively slow contrac¬ 
tion of the inner grid boundary. It would be highly desirable 
to perform model calculations also for faster boundary con¬ 
tractions, which is our plan for future work. In view of this 
caveat the arguments and test results discussed in this section 
should still be taken with a grain of salt. 

On grounds of the discussion of our results in the x-y-plane 
one can actually easily understand why the definition of the 
separation curve of exploding and non-exploding models in 


8 The radius R ns of the proto-neutron star is defined by the radial posi¬ 
tion where the density is 10 11 gem -3 . As in panels d-f of Fig.^] the time¬ 
averaging for (^ns) 3 (Jo ms * s performed from the passage of the s = 4 interface 
through the shock until the onset of the explosion for successful models and 
from the infall of the s = 4 interface until 300 ms later for non-exploding 
models. We also employ a normalization factor of 70 km to recover (roughly) 
the same range of values for x as in the case of x. 
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the present paper did not require us to take into account a 
possible dependence of the accretion luminosity on the NS 
radius. Instead, we could safely ignore such a dependence 
when we coined our ansatz that L v oc Mm'(r) oc M 4 // 4 . There 
are two reasons for that. On the one hand, the separation line 
in the M 4 / 44-//4 plane is fairly flat. A variation of the NS ra¬ 
dius corresponds to a horizontal shift of the location of data 
points in the x-y plane. However, with the contraction behav¬ 
ior of the neutron stars obtained in our simulations, only for 
(relatively few) points in the very close vicinity of the sepa¬ 
ration curve such a horizontal shift is sufficiently big to po¬ 
tentially have an influence on whether the models are classi¬ 
fied as non-exploding or exploding. On the other hand, the 
models near the separation line typically blow up fairly late 
(t eX p > 0.6 s), and the NS masses lie in a rather narrow range 
between roughly 1.4 M 0 and 1 .6 M 0 for the gravitational mass 
(see Fig.[l2|. For such conditions the variation of the NS ra¬ 
dius is ofsecondary importance (cf. Fig. 3 in Muller & Janka 


2014| >. Fow-mass progenitors with less massive N3s] whose 


radius (a t the sam e time) can be somewhat larger (Fig. 3 in 
|Miiller & Janka|2014) >, however, are located towards the lower 
left corner of the x-y plane and therefore far away from the 
separation curve (compare Figs.[TT|and[T 2 |. An incorrect hor¬ 
izontal placement of these cases (due to the omission of the 
dependence of the neutrino luminosity on the NS radius) does 
not have any relevance for the classification of these models. 

3.7. Brief comparison to previous works 

Our result of a complex pattern of NS and BH forming 
cases as function of progenitor mass was previously found 


by Ugliano et al. 

(2012 

jcha & Thompson 

(2015) 


too, and was confirmed by Pe- 
Aside from differences in details 


ferent SN 1987A calibration models, the main differences 


of the results presented here compared to those of Ugliano 
et al. (| 2012 |) are lower explosion energies for progenitors with 


M < 13 M q (see discussion in Sect. |2.3.4|> and lower fallback 


masses as mentioned in Sect. [33] Based on simple arguments 
(which, however, cannot account for the complex dynamics 
of fallback), Pejcha & Thompson (|2015j> already expected (in 
particular for their parameterization (a)) that cases with signif¬ 
icant fallback —in the sense that the remnant masses are sig¬ 
nificantly affected— should be rare for solar-metallicity pro¬ 
genitors. Our results confirm this expectation, although the 
ZAMS masses with significant fallback are different and less 


numerous than in the work by Pejcha & Thompson 
As pointed out by these authors, fallback has potentially - im¬ 
portant consequences for the remnant mass distribution, and 
the observed NS and BH masses seem to favor little fallback 
for the majority of SNe. 

As discussed in detail by Ugliano et al.| (2012 , our ex 


plosion models (as well as parameterization (a) otjPejcha & 
Thompson 2015) predicts many more BH formati on cases and 


more mass intervals of non-exploding stars than O'Connor & 
Ott| (|2011|>, who made the assumption that stars with com¬ 
pactness £ 2.5 > 0.45 do not explode. Roughly consistent 


with our results and those of 


Ugliano et al.]( j2012j i, [Horiuchi 


et al. (2014 1 concluded on the basis of observational argu 


ments (comparisons of the SN rate with the star formation 
rate; the red supergiant problem as a lack of Type-IIP SNe 
from progenitors above a mass of ~16M 0 ) that stars with an 
“average” compactness of 4 3.5 ~ 0-2 should fail to produce 
canonical SNe (cf. Fig. 15 in |Sukhbold et al.]2015| . 

A correlation of explosion energy and 5 ^Ni mass as found 


by Pejcha & Thompson] (j2015j) and [Nakamura et al.| (|2015 1 
(both, however, without rigorous determination of observable 
explosion energie s at infinity) and as suggested by observa¬ 
tions (see [Pejcha & Thompson||2015| for details) is also pre¬ 
dicted by our present results (but not by Ugliano et al.|2012 


with their erroneous estimate of the fallback masses and the 
overestim ated e xplosion energies of low-mass progenitors, 
see Sect. |2.3.4|i. Our models yield low explosion energies 
and low nickel production towards the low-mass side of the 
SN progenitors (see also Sukhbold et al. 20l3jt. In co ntrast, 
the modeling approach by Pejcha & Thompson (2015) pre¬ 
dicts a tendency of lower explosion energies and lower 56 Ni 
masses towards the high-mass side of the progenitor distribu¬ 
tion (because of the larger binding energies of more massive 
stars), although there is a large mass-to-mass scatter in all re¬ 
sults. This difference could be interesting for observational 
diagnostic s. The modeling approach by Pejcha & Thomp- 


(2015) seems to yield neutrino-driven winds that are con¬ 


siderably stronger, especially in cases of low-mass SN pro¬ 
genitors, than those obtained in the simulations of Ugliano 


et al. ( 2012 ) and, in particular, than those found in current, 
fully self-consistent SN and PNS cooling models, whose be¬ 
havior we attempt to reproduce better by the Crab-motivated 
recalibration of the low-mass explosions used in the present 
work. Moreover, for the 56 Ni-£ , e xD correlation reported by Pe¬ 
jcha & Thompson! (|2015|), a mass interval between ~14 M Q 
and ~15 M 0 with very low explosion energies and very low 
Ni production, whic h does not exist i n our models, also plays 
an important role. Nakamura et al. ( 2015| found a positive 
correlation of the explosion energy and the 56 Ni mass with 
compactness 43 . 5 - We can confirm this result for the nickel 


production but not for the explosion energy (cf. Sukhbold et 
|al.|2015[ Fig. 15 there). A possible reason for this discr epancy 


could be the consideration of “diagnostic” energies by Naka¬ 
mura et al. (2015} at model dependent final times of their sim- 


ulations - Instead of asymptotic explosion energies at infinity, 
whose determination requires seconds of calculation and the 
inclusion of the binding energy of the outer progenitor shells 
(see Fig. |3] and |Sukhbold et al.|20T5| , Fig. 6 , for the evolution 
of the energies in some exemplary Emulations). 

Since a comprehensive discussion of explosion energies, 
nickel production, and remnant masses is not in the focus of 
the present work, we refer th e reader interested in these as¬ 
pects to the follow-up paper by Sukhbold et al. (2015). For the 
same reason we also refrain here from more extended com¬ 
parisons to the progenitor-dependent explosion and remnant 


predic t ions o f other studies, in part icular th ose of O’Connor 
& Ott ( 2011 [ 1 , Nakamura et al.| < 2015), and |Pejcha & Th omp¬ 


son! ( |2015} . A detailed assessment of the different model- 
ing methodologies and their underlying assumptions would 
be needed to understand and judge the reasons for differences 
of the results and to value their meaning in the context of su¬ 
pernova predictions based on the neutrino-driven mechanism. 
Such a goal reaches far beyond the scope of our paper. 

4. CONCLUSIONS 

We performed ID simulations of SNe for a large set of 
621 progenitors of different masses and metallicities, includ¬ 
ing the solar-metallicity s2002 serie s (fWoosley et al.||200Z} 
previously investigated by Ugliano et aH ( 2012 > and the new 
s2014 models of Sukhbold & Woosley (2014) with their fine 
gridding of 0.1 M 0 in the ZAMS mass. 

In order to obtain SN explosions in spherical symmetry, we 
adopt the methodology of Ugliano et al.](|2012|), using ID hy- 
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drodynamics and approximate neutrino transport and a PNS- 
core neutrino source, but with improvements in a number of 
modeling aspects (e.g., a nuclear high-density EoS and a fully 
self-consistent implemen tatio n of nuclear burning through a 
small network, cf. Sect. 12. 1 b. Explosions are triggered by 
a neutrino luminosity that is sufficiently large to overcome 
the critical threshold condition for shock runaway. This lu¬ 
minosity is fed by a progenitor-dependent accretion compo¬ 
nent during the post-bounce shock-stagnation phase as well 
as a component from the high-density core of the nascent NS. 
The core emission is also progenitor-dependent, because it 
scales with the mass of the hot accretion mantle that assem¬ 
bles around the cooler high-density core of the PNS during 
the pre-explosion evolution. The conditions for an explosion 
are thus tightly coupled to the progenitor structure, which de¬ 
termines the post-bounce accretion history. 

Our approach contains a number of free parameters, whose 
values are calibrated by reproducing observational properties 
(explosion energy, 56 Ni mass, total neutrino energy and sig¬ 
nal duration) of SN 1987A with suitable progenitor models 
of this SN. We consider five different such progenitors for 
our study, namely besides the sl9.8 star of the s2002 model 
series of Woosley et al. ( 2002), which was used by Ugliano 


|et al.| l |2012) as calibration model, also 15, 18, and 20 M Q of 
Woosley and collaborators as well as a 20 M a model from 
Shigeyama & Nomotoj(jT990) (see Sect. |2.2|». 

Because ID simulations cannot properly reproduce the pe¬ 
riod of continued accretion and simultaneous outflow that 
characterizes the early expansion of the revived shock in 
multi-dimensional simulations and delivers the explosion en¬ 
ergy, our ID models exhibit an extended episode of accretion 
that is followed by a strong early wind phase. The overesti¬ 
mated mass loss during the latter phase compensates for the 
enhanced preceding accretion, and the associated recombina¬ 
tion energy yields the dominant contribution to the power sup¬ 
ply of the beginning explosion. A detailed discussion of our 


methodology can be found in Sect. 2.3 
Overall, our results confirm the ZAMS -mass dependent ex¬ 
plosion behavior that was found by Ugliano et al.|(j2012j) for 
the s2002 model series. For the same explosion calibration 
this progenitor set and the newer s2014 models have basic fea¬ 
tures in common. Moreover, for all calibration cases we ob¬ 
serve similar irregular patterns of successful explosions alter¬ 
nating with BH formation events above ~15 M Q . The largest 
fraction of BH formation cases is obtained with the w20.0 


calibration model, a 20 M 0 SN 1987A progenitor of Woosley 
et al. (l997j>, which explodes relatively easily and reproduces 


the SN 1987A 56 Ni yield with a fairly low ratio of explosion 
energy to ejecta mass of E exp /Mf,j ~ 0.7 only. The core neu¬ 
trino source is correspondingly weak and enables successful 
SNe in a smaller subset of progenitors. On the other hand, 
we obtain the closest similarity of the explodability of the 
investigated progenitor sets when we use the si9.8 (Ugliano 
et al.’s) calibration model and Shigeyama & Nomoto 1990 1 
n20.0 SN 1987A progenitor, which possess very similar com¬ 
pactness values ^ 2 . 5 - 

Non-exploding cases tend to correlate with local maxima 
of the compactness £ 2 . 5 , of the total binding energy outside 
of the iron core, E\,{m > Mp e ), and, most significantly, local 
maxima of the total binding energy E\,(m/M Q > M 4 ) outside 
of the mass coordinate M 4 = m(s = 4)/M G , where the dimen¬ 
sionless entropy per nucleon reaches a value of 4. Many (but 
not all) non-exploding progenitors below ~22 M Q also coin¬ 


cide with local minima of Mp e and, in particular, with minima 
of M4. However, there are no fix threshold values of any of 
these characteristic parameters of the pre-collapse progenitor 
structure that could be used to discriminate favorable from 
non-favorable conditions for an explosion. 

Guided by such insights we propose a two-parameter cri¬ 
terion to classify the explodability of progenitor stars by 
the neutrino-heating mechanism in dependence of the pre¬ 
collapse properties of these stars. The two structural param¬ 
eters that turn out to yield the best separation of successful 
and unsuccessful cases are M 4 and the mass derivative fj 4 = 
dm/dr[M G /(1000km)]~ 1 | i=4 = m'(s = 4)[M G /(1000km)r 1 
just outside of the s = 4 location, which we combine to a pa¬ 
rameter x = /V/4/./4. The parameters x and y = // 4 are tightly 
connected to the two crucial quantities that govern the physics 
of the neutrino-driven mechanism, namely the mass-accretion 
rate of the stalled shock, M, and the neutrino luminosity L v . 
The former determines the ram pressure that damps shock ex¬ 
pansion and can be mathematically linked to the mass deriva¬ 
tive m' (see Eq. [5]). The latter is a crucial ingredient for the 
neutrino heating that is responsible for shock revival and is 
dominated by the accretion luminosity and the PNS-mantle 
cooling emission during the crucial phase of shock revival. 
Both of these scale with M and/or the accretor mass (i.e., the 
PNS mass) so that L v oc M ns M expresses the leading depen¬ 
dence. Since the neutrino-driven explosions in our simula¬ 
tions set in shortly after the entropy interface and density jump 
around s = 4 have fallen through the shock (Fig. [2]), M4 can 
be taken as a good proxy of M ns as the accretor mass, and /./ 4 
can serve as a measure for the mass accretion rate M of the 
PNS. 

Higher values of M4 tend to be favorable for an explosion 
as shown by Fig. [9] where many non-exploding cases (dark 
gray) correlate with local minima of M 4 . The reason is that 
the neutrino luminosity scales (roughly) with x = /W 4/./4 (the 
actual sensitivity of the neutrino-energy deposition to M4 is 
even steeper). Therefore higher M4 imply greater neutrino lu¬ 
minosities and stronger neutrino heating. In contrast, the in¬ 
fluence of y = jj 4 is ambivalent. On the one hand a high value 
of /j 4 increases the neutrino luminosity, on the other hand it 
also causes a large ram pressure that has to be overcome by 
neutrino heating. The effect of these competing influences is 
that a higher value of M4 in association with a lower value 
of fi4 favors explosions. Reversely, non-exploding cases in 
Fig. 0 are correlated with local minima of M4 and maxima 
of H4. Visually, this means that explosion cases are preferen¬ 
tially located toward the lower right of the M4H4-/44 parameter 
space (cf. Figs.|6][8|[T2). 

The parameters x and y therefore span a plane in which suc¬ 
cessful explosions and failures with BH formation are clearly 
separated. The progenitor stars populate an astonishingly nar¬ 
row band that stretches from the lower left corner with the 
lightest stars to the upper right direction of the x-v-plane, 
where the massive progenitors with the biggest iron cores and 
highest binding energies of overlying material are located (see 
Fig. [IT}. While SNe can be found in the region of low values 
of y, i.e., for low mass-accretion rate, the non-exploding cases 
inhabit the domain of high y, but the limiting value of the 
mass-accretion rate that prevents the success of the neutrino- 
driven mechanism grows with the value of x. Both sectors in 
the x-y plane are separated by a boundary line that can be well 
represented by a linear function y sep (x) (Eq. [9ji with increas¬ 
ing slope. (The values of the dimensionless coefficients of 
this linear relation are listed for all calibration models in Ta- 
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ble[2]) Because of the physical meaning of the parameters x 
and y, there is a close correspondence between the separation 
line y sep (x) and the critical threshold luminosity L v ,crh(M) that 
has to be exceeded to trigger runaway expansion of the accre¬ 
tion shock by neutrino heating. The rising slope of the sep¬ 
aration curve in this context means that for each value of the 
neutrino luminosity, respectively x, there is an upper limit of 
the mass accretion rate, respectively y, up to which neutrino- 
driven explosions are possible, and that this BH-formation 
threshold value of y grows with x. The parameters x and y, 
computed from the progenitor profiles before collapse, allow 
one to judge whether a considered star is able to overcome 
the threshold neutrino luminosity for an explosion, or, in other 
words, whether its mass accretion rate stays below the critical 
limit above which the onset of the explosion is prevented. 

Only ~l-3% of all investigated progenitors do not follow 
this discrimination scheme in their final fate but lie on the 
wrong side of the separation curve. These outliers are char¬ 
acterized by pathologies of their entropy and density pro¬ 
files that describe the composition-shell structure in the Si-O- 
layers. Such special conditions lead to unusual combinations 
of mass accretion rate and PNS masses. Our two-parameter 
criterion expressed by the separation line y cn i(x) therefore en¬ 
ables one, with a very high significance, to predict the explod- 
ability of progenitor stars via the neutrino-driven mechanism 
by referring to the properties of these stars as captured by the 
pair of parameters x and y. 

Clausen et al. ( j2015j ) explore the interesting possibility 
that the death of massive stars in NS vs. BH formation may 
be better captured by a probabilistic description. The non¬ 
monotonic variations of explosion vs. non-explosion with 
ZAMS mass or compactness might be interpreted as stochas- 
ticity in the explosion behavior. However, by considering the 
problem in a more appropriate two-parameter space, our two- 
parameter criterion unmasks these putatively random varia¬ 
tions as actually deterministic phenomenon. Clausen et al. 
(2015), in contrast, suggest a variety of factors besides ZAMS 


mass and metallicity, e.g. rotation, binarity, the strength of 
magnetic fields, stochastic differences in the pre-collapse 
structure or even in the explosion mechanism, that might in¬ 
troduce a randomness such that a star of given mass might not 
form either a NS or a BH but both with a certain probability. If 
such a diversity in the stellar destiny depends on a causal pro¬ 
cess, for example the presence of different amounts of spin, 
a deterministic description could still apply but would require 
an extension to a parameter space of more than the two dimen¬ 
sions combined by our current criterion, e.g., by adding extra 
dimension(s) that capture the role of rotation in the explosion 
mechanism. If, in contrast, truly stochastic effects like turbu¬ 
lent processes or chaotic fluctuations in the progenitor, decide 
about the stellar fate, a deterministic criterion for explodabil- 
ity would be ruled out and a probabilistic description would 
be indispensable. 

Our study has a number of more caveats that need to 
be addressed. The understanding of the explosion mecha- 
nism(s) of massive stars is still incomplete, although con¬ 
siderable progress has been achieved in recent years due to 
the progress in 2D and 3D modeling and in particular also 
through improvements in the treatment of the crucial neutrino 


lanka et al. 2012} lanka 

2012; Burrows 2013 Mezzacappa 

et al. 2015 ; Melson et al. 

2015afb; Eentz et al. 2015 and ref- 


erences therein). Without self-consistent 3D explosion sim¬ 
ulations for large sets of progenitor stars being possible yet, 


our study refers to a ID modeling approach, in which not only 
the neutrino description is approximated in many aspects, but 
also the explosions have to be triggered artificially. We em¬ 
ploy a boundary condition that replaces the high-density, low- 
entropy core of the nascent NS as a neutrino source, and we 
describe the time-dependent behavior of this core and of the 
coupling between the core and its hot accretion mantle by a 
simple, analytical model. Calibrating the involved free pa¬ 
rameters by observational constraints from SN 1987A for the 
more massive stars and by comparison to results of sophis¬ 
ticated SN models for low-mass progenitors is intended to 
anchor our simplified description on empirical ground. Al¬ 
though the elements of this approximate approach appear 
qualified to capture the essence of the neutrino-heating mech¬ 
anism in dependence of the prog enitor-specific post-bounce 
accretion evolution (cf. Sect. |2.3| for details), confirmation by 
fully self-consistent, multi-dimensional SN calculations is ul¬ 
timately indispensable. It is also evident that our study, which 
is only concerned with neutrino-driven explosions, cannot 
yield any information about the possibility and implications 
of other mechanisms to blow up stars, for example magne- 
torotational explosions of rapidly rotating stellar cores, which 
might be a consequence of magnetar or BH formation and 
could be associated with hypernovae and gamma-ray burst 


SNe (see, e.g., Mazzali et al. 2014 

i as well as ultra-luminous 

SNe (WoosleypOlO; Kasen & Bile 

sten 2010; Sukhbold et al. 


[2QI5F 


Our study employs pre-collapse models that emerge from 
ID stellar evolution calculations of single, non-rotating SN 
progenitors with prescribed mass-loss histories. The results 
of our study naturally depend on the post-bounce accretion 
properties of the collapsing stars. The (iron and low-entropy) 
core masses as well as the entropy and density jumps at the 
composition-shell interfaces play an important role in setting 
the conditions for the neutrino-heating mechanism, which is 
obvious from the definition of our parameters x and y. It 
is conceivable that multi-dimensional hydrodynamics could 
lead to considerable changes of the stellar propertie s as func¬ 
tions of the progenitor mass (e.g., Arnett et al. 2015| ), and 
that asymmetries and perturbations in the shell-burning layers 
of the pre-collapse core might have important consequences 
for the development of SN explosions by the neutrino-driven 
mechanism (e.g., Arnett & Meakin|201 lUCouch & Qtt|2013 
... ’ ' " ' ' We 


|2015| [Muller & Janka||2015[ |Couch et al.||2015| l. We are 
hopeful that the basic insights of our study, in particular the 
existence of a two-parameter criterion for the explodability 
—expressed by a SN-BH separation line y cr itM in the x-y- 
space and the tight connection between this curve and the 
critical luminosity condition L vcnl (M) of the neutrino-driven 
mechanism— possess more general validity. The explosion 
properties of the progenitor stars as functions of the ZAMS 
mass, however, do not only depend (moderately) on the con¬ 
sidered SN 1987A progenitor models but will probably also 
change once multi-dimensional stellar evolution effects will 
be accounted for in the pre-supernova conditions. 
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